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ABSTRACT 


The  correlation  of  meteorological  events  such  as  the  jet  stream,  gravity  waves  and 
boundary  layer  circulation  with  the  optical  turbulence  parameters,  the  transverse  co- 

*  herence  length  r0  and  the  isoplanatic  angle  0O  is  essential  for  interpreting  and  forecasting 
imaging  and  laser  systems  performance.  In  support  of  the  United  State  Air  Force  Relay 
Mirror  Experiment,  the  Naval  Postgraduate  School  performed  a  series  of  six  site  char¬ 
acterization  measurements  near  Kihei,  Maui  during  August  1987  to  July  1988.  Spatial 
and  temporal  summaries  of  atmospheric  events  corresponding  to  the  optical  remote 
sensor  data  are  presented  using  meteorological  data  from  the  National  Weather  Service 
Rawinsonde  Observation  Stations,  synoptic  charts,  GOES- WEST  infrared  satellite  im¬ 
ages  and  four  Kihei,  Maui  rawinsonde  datasets.  To  quantify  the  correlation  between 
optical  turbulence  measurements  and  meteorological  phenomena,  four  methods  of  cal¬ 
culating  Cf  from  rawinsonde  data  were  investigated.  Results  show  that  existing 
rawinsonde  systems  are  inadequate  for  direct  Or  calculation.  However,  moderate  im- 

•  provements  in  the  vertical  resolution,  the  temperature  resolution  and  probe  response 
time  will  allow  direct  calculations  of  optical  turbulence  parameters  from  rawinsonde 
data. 
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that  at  least  some  of  the  many  blessings  I've  received  while  pursuing  the  thesis  have  also 
been  enlightening  and  encouraging  to  those  with  whom  I  have  run  this  race. 


I.  INTRODUCTION 


Atmospheric  turbulence  produces  optical  phase  and  amplitude  fluctuations  that  se¬ 
riously  degrade  laser  beam  propagation  and  the  performance  of  imaging  systems.  In 
order  to  perform  systems  analyses  of  optical  propagation,  the  spatial  and  temporal  dis¬ 
tribution  of  atmospheric  optical  turbulence  must  be  understood.  Over  the  past  two  and 
a  half  years  optical  turbulence  measurements  in  the  form  of  isoplanatic  angles,  0O,  and 
transverse  coherence  lengths,  r0,  have  been  made  at  a  variety  of  sites  in  the  continental 
United  States,  as  well  as  Maui,  Hawaii  and  Vieques,  Puerto  Rico.  Though  a  dominant 
application  of  these  measurements  involve  laser  technology  and  astronomy,  the  inter¬ 
pretation  of  the  results  are  very  valuable  to  the  meteorologist.  The  optical  data  repres¬ 
ents  an  instantaneous  summary  of  the  optical  turbulence  along  a  vertical  path. 
Imbedded  within  this  vertical  cross-section  is  information  describing  various  synoptic 
and  mesoscale  phenomena  such  as  the  jet  stream,  gravity  waves,  low  level  jets  and 
boundary  layer  circulation.  By  correlating  meteorological  measurements  and  optical 
turbulence  parameters,  the  practical  cause  and  effect  relationships  can  be  assimilated 
into  prognostic  tools. 

The  United  States  Air  Force  Relay  Mirror  Experiment  (RME)  requires  a  critical 
evaluation  of  atmospheric  turbulence  over  the  Maui,  Hawaii  experiment  site.  RME  is 
designed  to  "assess  how  accurately  an  orbiting  mirror  can  place  (a)  reflected  (laser)  beam 
on  a  ground  target  with  low  beam  wander"  (Nelson,  1988).  Observed  beam  jitter  and 
wander  arc  a  product  of  both  the  laser  system  and  the  atmospheric  turbulence.  Accurate 
measurements  of  the  atmospheric  contributions  are  vital  for  evaluating  the  optical  sys¬ 
tem's  performance. 

Between  August  1987  and  July  1988,  the  Naval  Postgraduate  School  (NPS)  per¬ 
formed  six  measurement  sessions  of  approximately  one  week  duration  each  in  order  to 
quantify  the  RME  target  board  site's  optical  turbulence  characteristics.  NPS  optical 
instruments  were  used  to  acquire  atmospheric  coherence  length,  r0,  and  isoplanatic  angle, 
0O,  measurements  near  sea  level  at  Kihei,  Maui.  In  addition  to  meteorological  data  from 
the  National  Weather  Service  Rawinsonde  Observation  Stations,  synoptic  charts  and  the 
GOES-WEST  Infrared  Satellite  Images,  four  of  the  six  Kihei  data  sets  had  on-site 
rawinsonde  measurements. 
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This  thesis  attempts  to  correlate  the  RME  meteorological  measurements  with  the 
optical  data.  Initially,  the  National  Weather  Service  measurements,  charts  and  satellite 
images  were  used  to  compare  meso-  and  synoptic  scale  events  with  the  optical  results. 
This  method  provided  qualitative  results. 

The  major  thrust  of  this  thesis,  however,  materialized  with  the  availability  of  high 
resolution  rawinsonde  soundings.  These  profiles  finally  provided  the  resolution  neces¬ 
sary  for  attempting  the  direct  computation  of  optical  parameters  from  the  rawinsonde 
data.  The  August  and  December  1987  sessions  used  the  VIZ  WO-8000RP+  rawinsonde 
system.  With  a  vertical  resolution  of  approximately  35  m,  measurements  were  recorded 
up  to  12  km.  The  April  and  July  1988  sessions  made  use  of  the  Vaisala  DigiCORA  MW 
1 1  system.  This  system  profiled  meteorological  data  up  to  approximately  24  km  with  a 
sampling  rate  of  2  seconds  (surface  to  2  km)  and  10  seconds  (above  2  km).  The  8  m  and 
40  m  vertical  resolutions  of  the  DigiCORA  system  were  considered  reasonable  for  com¬ 
puting  the  optical  turbulence  directly  from  the  rawinsonde  data. 

The  optical  variables  measured,  r0  and  0O,  depend  on  an  integral  of  the  index  of  re¬ 
fraction  structure  parameter,  Q,  along  the  vertical  path.  The  parameter  Q  is  propor¬ 
tional  to  the  temperature  structure  parameter  Ot  which  can  be  computed  from  the 
rawinsonde  data.  Four  methods  for  computing  the  temperature  structure  parameter  Or 
from  the  rawinsonde  data  were  investigated: 

1.  the  Simple  Potential  Temperature  Differences  (SPTD)  Equation, 

2.  the  Potential  Temperature  Differences  with  a  Running  Mean, 

3.  the  Bulk  Gradient  Equation  (BG)  and 

4.  conditional  data  sampling  based  on  the  Richardson  Number. 

Each  Or  method  was  tested  by  calculating  r0  and  0O  from  the  Ot  profile.  These  results 
were  ultimately  compared  to  the  optical  measurements  simultaneously  acquired  over  the 
same  site. 
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II.  BACKGROUND 


A.  HISTORICAL  REVIEW  OF  CORRELATION  RESEARCH 

International  scientific  interest  in  the  meteorological  influence  on  "seeing",  a  rating 
system  based  on  the  angular  size  of  a  star  image,  flourished  in  the  mid  to  late  50's. 
Following  a  brief  decline,  interest  was  reestablished  in  1962.  The  new  upper  atmospheric 
and  outer  space  environmental  data  introduced  by  artificial  satellites  and  space  probes 
provided  a  catalyst  for  the  development  of  atmospheric  physics.  One  of  the  dominant 
questions  addressed  by  the  melange  of  scientific  backgrounds  required  in  the  atmo¬ 
spheric  physics  field  involved  how  meteorological  conditions  affect  "seeing"  (Kucherov, 
1966a). 

In  October  1962,  the  issue  of  observatory  site  selection  and  related  astroclimatic 
topics  were  discussed  at  the  International  Astronomical  Union  Symposium  held  in 
Rome.  Geophysicists,  meteorologists  and  astronomers  worked  in  cooperation  to  study 
the  interrelationships  between  atmospheric  conditions  and  astronomical  "seeing".  Par¬ 
ticipants  resolved  that  placing  telescopes  at  certain  heights  above  the  ground  would  im¬ 
prove  image  quality  (Kucherov,  1966a).  The  following  summarizes  some  of  the 
pertinent  articles  presented  at  the  Third  All-Union  Conference  held  four  months  earlier 
in  Kiev. 

Kucherov  (1966b)  observed  "shimmer  angle"  and  meteorological  conditions  at  one 
seashore  and  two  inland  sites  over  a  period  of  approximately  two  years  (1960-1962). 
The  results  of  that  study  confirmed  the  relationship  between  synoptic  fronts  (mostly  cold 
fronts)  and  astronomical  "seeing".  It  should  be  noted  that  25%  of  the  cases  display  no 
relationship. 

Bol'shakova  and  Darchiya  (1966)  investigated  the  causes  of  "seeing"  condition 
changes.  While  they  acknowledge  the  synoptic  fronts  as  an  influence  on  "seeing",  they 
suggest  that  "seeing"  is  affected  by  frontal  passages  only  when  air  masses  with  optical 
characteristics  differing  from  the  original  advect  over  the  site.  Using  statistics,  they 
support  their  claim.  To  a  first  approximation,  50%  of  the  cases  show  the  synoptic  front 
causing  a  deterioration  in  "seeing",  or  alternatively,  absence  of  a  synoptic  front  ensures 
high  "seeing"  conditions.  The  rest  of  the  cases  display  the  reverse  situation.  In  their 
opinion,  the  analyois  must  be  repeated  taking  into  account  the  type  of  front,  the  incli- 


nation  of  the  frontal  plane,  the  direction  of  the  front's  motion,  and  ideally,  other  detailed 
meteorological  information  defined  in  various  atmospheric  layers. 

Zinchenko  (1966)  comes  to  the  following  conclusions: 

1.  The  lower-atmospheric  turbulence  is  of  decisive  significance  in  all  that  concerns 
"seeing": 

a.  "optical  inhomogeneities  are  most  effective  in  the  ground  layer,  where  the  air 
density  is  highest;" 

b.  "turbulence  intensity  is  greatest  in  this  layer";  and 

c.  "eddies  of  minimum  size  occur  in  the  ground  layer". 

2.  Contradictory  results  exist  with  the  applicability  of  Richardson  number(Ri)  as  a 
criterion  of  high  atmospheric  turbulence. 

The  main  reason  supporting  the  latter  point  is  that  turbulence  intensity  is  not  a 
single-valued  function  of  Ri,  though  a  dependency  on  Ri  exists.  Zinchenko  explains, 
when  Ri  <  1,  one  can  qualitatively  accept  that  turbulence  set  in;  and,  the  amount  of  ac¬ 
cumulated  eddy  energy  in  steady-state  conditions  increases  with  the  decrease  in  Ri. 
(Note:  Current  understanding  requires  Ri<0.25  for  the  above  to  be  true.)  Though  lack 
of  data  prohibited  an  exact  definition  of  the  layer  directly  responsible  for  the  blurring 
of  the  stellar  diffraction  disks,  Zinchenko  found  that  turbulence  does  show  a  tendency 
to  be  restricted  to  the  0. 5-2.0  km  layer. 

B.  THE  RELAY  MIRROR  EXPERIMENT 

The  Relay  Mirror  Experiment  (R.Y1E)  is  a  United  States  Air  Force  Experiment 
monitored  by  the  Kirtland  Air  Force  Weapons  Laboratory,  Albuquerque,  New  Mexico. 
RME  is  designed  to  "assess  the  ability  to  place  a  laser  beam  on  a  ground  target  with  low 
beam  wander"  (Nelson,  198S).  As  jescribed  by  Nelson,  RME  involves  the  propagation 
of  three  ground  to  space  (uplink)  laser  beams  as  well  as  the  reflection  from  space  to 
ground  (downlink)  of  one  of  these  beams  (see  Figure  1).  Two  of  the  uplink  lasers,  one 
from  AMOS  (Air  Force  Maui  Observatory  Station),  a  site  located  at  the  Halcakala 
crater  rim,  and  the  other  from  Kihei,  a  site  located  at  the  volcano's  base,  serve  to  align 
a  relay  mirror  located  on  a  payload  orbiting  over  the  Maui,  Hawaii  stations.  With  the 
mirror  aligned,  a  1.06  fx  m  laser  travels  from  AMOS  to  payload  and  is  relayed  down  to 
a  hexagonal  target  board  at  Kihei.  The  telescopes  and  optical  detectors  on  the  target 
board  record  the  centroid  motion  and  beam  size  of  the  reflected  1.06  n  m  laser  beam. 
Atmospheric  turbulence  has  an  important  effect  on  both  uplink  and  downlink  beams 
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(Nelson,  1988).  Understanding  and  characterizing  the  atmospheric  turbulence  over  the 
Hawaiian  experiment  site  is  crucial  for  assessing  the  success  of  the  RME  project. 

C.  HAWAIIAN  CLIMATOLOGY 

The  following  is  a  brief'summary  of  Hawaiian  climatology  with  a  focus  on  the  main 
experiment  site  in  Kihei,  Maui. 

The  tropical  island  of  Maui  is  located  at  156°  15'  longitude  and  20°  45'  latitude 
(Blumenstock  and  Price,  1967).  The  island  is  the  second  most  southern  island  in  the 
Hawaiian  chain  and  has  an  extreme  length  of  77  km,  maximum  width  perpendicular  to 
extreme  length  of  42  km  and  a  maximum  elevation  of  3. 1  km.  Nicknamed  "The  Valley 
Isle",  this  1886  square  km  island  is  comprised  of  an  isthmus  linking  two  volcanos.  The 
older  volcano,  Pui  Kukui,  is  in  the  northwest  sector  reaching  a  peak  of  1.8  km.  To  the 
southeast  is  the  dormant  volcano,  Haleakala.  Though  Halcakala  defines  the  island's 
highest  point,  it's  crater  is  one  of  the  largest  on  earth:  12  km  long,  4  km  wide  and  about 
900  m  deep  (Creal  et  al.,  1987).  Topographically  induced  high  winds  crossing  the  sum¬ 
mit  arc  common. 

1.  Climatic  Regions  of  Hawaii 

The  tropical  climate  on  Maui  varies  greatly  with  location.  The  dominate  sub¬ 
climates  include  the  following: 

1.  Ralnv  MOUNTAIN  slopes  on  the  north,  northeast  WINDWARD  SIDE  (above 
600  m): 

2.  Moderately  rainy  and  mild  WINDWARD  LOWLANDS  (below  600  m)  to  the 
north,  northeast; 

3.  Dry  LEEWARD  LOWLANDS  to  the  south,  southeast; 

4.  Sometimes  rainy,  sometimes  dry  INTERIOR  LOWLANDS  between  the  two 
volcanic  peaks; 

5.  Dry,  cold  HIGH  MOUNTAINS  above  600-900  m  (Blumenstock  and  Price,  1967). 

Rainfall  measurements  reveal  the  sharp  weather  contrasts  over  this  island.  On 
the  leeward  coastal  areas  such  as  Kihei,  and  mountain  summits,  average  annual  rainfall 
is  0.5  m.  The  lower  windward  slopes  of  a  high  mountain  and  at/near  the  summit  of 
lower  mountains  (western  Maui)  have  an  average  annual  rainfall  of  7.6  m  (Blumenstock 
and  Price,  1967). 

2.  Circulation 

The  dominant  atmospheric  circulation  for  the  Hawiiain  islands  is  the  north- 
northeasterly  trade  winds.  The  prevalence  of  this  flow  reflects  the  seasonal  north-south 
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EXPERIMENT  SYSTEM  SCHEMATIC 
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.  1.  Relay  Mirror  Experiment  System  Schematic 
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movement  of  the  North  Pacific  High.  From  May  to  September  (summer),  the  North 
Pacific  High  is  to  the  north.  During  this  season,  the  southern  sectors  of  the  counter¬ 
clockwise  flow  around  the  High  is  such  that  trades  prevail  80-95%  of  the  time.  From 
3  October  to  April  (winter)  the  North  Pacific  High  drops  south  such  that  the  heart  of  the 

trades  pass  well  to  the  south  of  the  islands.  The  frequency  of  trades  during  this  season 
drops  to  50-S0%  (Blumenstock  and  Price,  1967). 

Land  and  sea  breeze  circulation  occurs  on  a  much  smaller  scale.  Kihei  is  unique 
in  that  land  and  sea  breezes  are  present  even  during  moderate  to  strong  trades 
(Blumenstock  and  Price,  1967).  The  sometimes  brisk  sea  breezes  occur  from  forenoon 
to  early  evening.  During  the  night,  a  very  gentle  land  breeze  flows  until  shortly  after 
sunrise.  The  Kihei  region  also  receives  a  nighttime  Katabatic  (downslope)  flow  compli¬ 
menting  the  land  breeze  already  in  progress  (personal  observation). 

3.  Temperature,  Humidity,  Cloud  Cover 

Hawaii's  mean  annual  temperature  range  generally  fluctuates  less  that  6  degrees 
Celsius  (Blumenstock  and  Price,  1967).  Daily  summer  temperatures  at  2.5  km  are  be¬ 
tween  4C  (nighttime)  and  10C  (daytime).  The  winter  range  at  this  height  is  0-10C.  In 
the  lowlands,  daytime  temperature  is  between  21-27C,  while  nighttime  is  in  the  16-2 1C 
range.  The  warmest  months  in  the  dry  areas  such  as  Kihei  are  August  and  September 
(Blumenstock  and  Price,  1967). 

Generally  humid  conditions  exist  with  high  clouds  except  in  the  driest  coastal 
regions  and  at  high  elevations.  In  the  absence  of  synoptic  scale  events,  Kihei  begins  a 
day  with  locally  clear  skys.  A  backdrop  of  clouds  clinging  to  the  north  and  west  slopes 
of  Haleakala  is  common.  Throughout  the  day  Kihei  experiences  a  growing  and  decaying 
of  cloud  cover.  From  sunset  to  approximately  2000  hours  local  time,  the  partly  cloudy 
sky  cover  over  the  site  dissipates.  Because  of  the  dramatic  climate  changes  over  such  a 
small  island,  clouds  arc  generally  visible  at  all  hours  from  the  Kihei  experiment  site  in 
some  section  of  the  sky  (Blumenstock  and  Price,  1967). 

4.  Storms 

Generally,  Hawaiian  weather  falls  into  trade  wind  and  non-trade  wind  condi¬ 
tions.  During  trade  wind  conditions  one  can  expect  trade  wind  showers  on  the  E  and 
NE  (windward)  coasts;  possible  light  trade  wind  showers  over  the  S  and  SW  (leeward) 
areas;  and,  sea  breezes  ascending  mountain  slopes  potentially  producing  light  to  occa- 

r 

sionally  heavy  showers  in  the  late  afternoon  and  evening.  Non-trade  wind  conditions 
are  associated  with  light  and  variable  winds;  moist  air  slowly  advecting  onto  the  islands 
from  the  southeast,  south  or  southwest;  intense  local  storms;  and,  the  possibility  of 
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major  and  minor  storm  systems  passing  through  the  islands  (Blumenstock  and  Price, 
1967). 

Major  storms  for  Hawaii  fall  into  one  of  four  classes  of  disturbances: 

1.  Cold  Front  Storms  -  winter  storm  systems  associated  with  cold  fronts  that  ap¬ 
proach  from  the  northwest  end  of  Hawaii,  and  move  eastward  barely  reaching 
Oahu. 

2.  Kona  Storms  -  winter  storm  systems  similar  to  cold  front  storms;  winds  are  gener¬ 
ally  from  "kona"  or  leew'ard  directions  relative  to  the  trades;  well-developed  kona 
storms  are  more  widespread  and  more  prolonged  than  the  usual  cold  front  storms. 

3.  Tropical  Storms  -  similar  to  hurricanes  but  with  winds  below  120  km/hr. 

4.  Hurricane  -  the  least  prevalent  storm  system;  in  63  years  only  four  hurricanes  have 
crossed  over  the  islands.  Unlike  the  Cold  Front  and  Kona  Storms,  the  Tropical 
Storm  and  Hurricane  can  occur  at  any  tame.  They  are  most  likely  to  occur,  howr- 
ever,  between  the  months  of  July  to  December  (Blumenstock  and  Price,  1967). 

These  events  will  be  discussed  in  greater  detail  as  they  pertain  to  a  specific  case  study. 
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III.  ATMOSPHERIC  OPTICS  THEORY 


A.  ATMOSPHERIC  INFLUENCE  ON  LIGHT  PROPAGATION 

Light  propagates  through  a  medium  in  the  form  of  a  wavefront,  "a  surface  over 
which  an  optical  disturbance  has  a  constant  phase"  (Hecht  and  Zajac,  1974).  Fermat's 
Principle  explains  that  the  route  light  travels  will  be  the  smallest  Optical  Path  Length 
(OPL).  This  is  expressed  mathematically  as  (Hecht  and  Zajac,  1974) 

OPL=  [  n(s)ds ,  (1) 

where  n  is  the  index  of  refraction  (the  "speed  of  an  electromagnetic  wave  in  vacuum  to 
that  in  matter"),  and  s  is  the  path  length.  The  phase  of  a  propagating  wave  will  be 
k  —  lnl).  times  OPL  (Eq.  (1)).  Optical  turbulence  arises  when  the  accumulated  phase 
along  adjacent  paths  is  different.  Phase  distortion  in  the  atmosphere  occurs  when  a 
wavefront  encounters  random  irregularities  in  the  atmosphere's  index  of  refraction. 
Such  irregularities  are  due  to  the  turbulent  mixing  processes  of  the  atmosphere.  Thus, 
when  a  wavefront  from  a  laser  propagates  vertically  through  the  atmosphere,  accumu¬ 
lated,  random  phase  differences  will  degrade  the  laser  beam  and  image  system  perform¬ 
ance.  Typical  atmospheric  distortions  result  in:  laser  beam  centroid  wander, 
scintillation,  image  breakup  and  blurring. 

The  study  of  these  atmospheric  random  turbulence  processes  requires  the  utilization 
of  statistical  quantities,  means,  variances  and  high  moments.  Frequently,  variances  and 
the  higher  order  moments  are  not  single-valued,  unique  or  even  bounded.  Also,  the  at¬ 
mosphere  is  neither  homogeneous  nor  isotropic.  Kolomogorov  and  Tatarski  introduced 
the  concept  of  structure  functions  and  local  homogeneity  to  resolve  these  problems 
(Tatarski,  1961). 

Kolmogorov  and  Tatarski  define  the  structure  function,  DR{r)  ,  for  an  arbitrary 
spatially  distributed  random  variable  R(r)  as 

£>,(?)  =  <(/!(?,)  -K(f2)]2>,  (2) 

where  r  =  f,  —  r2 ,  and  the  angled  brackets  indicate  an  ensemble  average  of  the  enclosed 
quantity.  Assuming  the  random  process  generating  R  is  homogeneous  and  isotropic, 
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then  only  the  ma£.r,:'  ude  of  r  is  important,  r  =  |F,  —  F,|.  Local  homogeneity  assumes  that 
the  mean  <tf(r,)>  =  <R(r2)>  (Clifford,  1978).  Restricting  the  separation  r  to  the  inertial 
subrange  of  turbulence,  Kolmogorov  (1961)  develops  a  universal  form  of  the  single 
structure  function 


DH(r)-Clr213, 


(3) 


where  Q  is  a  structure  parameter,  a  measure  of  "the  total  amount  of  energy  in  the  tur¬ 
bulence"  (Clifford,  1978);  r  is  the  distance  between  the  two  sampled  points;  and,  r  must 
be  between  the  turbulence  innerscale,  /„,  and  the  turbulence  outerscale,  .  Eddy  sizes 
below  the  innerscale  correspond  to  the  eddy's  energy  dissipation  through  viscous  effects. 
The  outerscale  corresponds  to  "the  largest  scale  size  for  which  the  eddies  may  be  con¬ 
sidered  to  be  isotropic”  (Clifford,  1978). 

For  optical  propagation  through  turbulence,  the  characterizing  quantity  is  the  index 
of  refraction  structure  parameter  (Walters  and  Kunkel,  1981) 


,2  <(/?)  ~  /?;)  > 

“  ..2/3 


(4) 


where  //,  =  «(/-,)  and  n2  —  ;t(r2)  arc  the  atmospheric  index  of  refraction  at  points  r,  and  rh 
respectively.  Though  Q  is  difficult  to  measure  directly,  profiles  calculated  from 
rawinsondc  sensor  measurements  show  that  Q  generally  decreases  with  altitude  (Nelson, 
198S).  A  more  practical-approach  for  measuring  atmospheric  turbulence  is  to  measure 
the  temperature  structure  parameter  C|, 


.2  <T ,  -  T2f> 

C' - ^3  • 

Since  the  atmospheric  index  of  refraction 

n-  1  =  (79jd0'6)y 


(5) 


(6) 


depends  on  pressure  P  and  temperature  T  (Tatarski,  1961),  Q  is  related  to  Or  by  sub¬ 
stituting  Eq.  (6)  in  Eq.  (4)  to  get 


C/j  =  (79  x  10-6  )2Cf, 


(7) 
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assuming  turbulence  is  isobaric.  Specific  methods  for  calculating  Ot  will  be  discussed  in 
Chapter  5. 

B.  OPTICAL  VARIABLES 

As  described  earlier,  the  laser  beam  is  affected  by  atmospheric  turbulence  as  a  result 
of  index  of  refraction  fluctuations  along  the  turbulent  propagation  path.  Two  parame¬ 
ters  that  characterize  the  optical  turbulence  generated  by  the  atmosphere  are:  the 
transverse  coherence  length,  /•„,  and  the  isoplanatic  angle,  60 . 

1.  Transverse  Coherence  Length 

The  transverse  coherence  length,  r0,  is  a  measure  of  the  lateral  autocorrelation 
length  of  the  electric  field  propagating  through  the  turbulent  atmosphere.  Because  of 
the  Fourier  transform  properties  of  an  imaging  system,  r0  is  also  a  measure  of  the  spatial 
frequency  response  of  the  atmosphere,  the  atmospheric  modulation  transfer  function 
(MTF).  Viewed  as  a  component  of  a  linear  system,  the  atmosphere  acts  like  a  low  pass 
filter  suppressing  the  high  spatial  frequencies,  or  details  of  an  image. 

The  perturbation  theory  developed  by  Tatarski  and  extended  by  Fried  (1966) 
and  Lutomirski  and  Yura  (1971)  shows  that  the  long  term,  modulation  transfer  function 
of  the  atmosphere  is  expressed  in  terms  of  the  wave  structure  function  D  by 

MTF(v)  =  exp(  ),  (8) 

4m 

where  v  is  the  spatial  frequency  in  cyclcs/m,  /  is  the  wavelength  of  light  and /is  the  focal 
length  of  the  optical  system.  For  a  plane  wave,  D  has  the  form 


Z)(p)  =  2.91/cV/3 


Cl 


Cfc)dzt 


Jo 


(9) 


where  k  is  2n\).,  the  wavenumber;  p  =  |F,  -  f2|  =  Xfv,  a  distance  measured  at  the  entrance 
pupil  of  the  optical  system,  L  is  the  optical  path  length  and  Q  is  the  index  of  refraction 
structure  parameter. 

The  general  format  of  MTF(v)  is  exp(  -ap!li)  where  a  is  a  turbulence  and  wave¬ 
length  scaling  parameter.  MTF{v)  is  frequently  written  as: 

MTF(r)  =  exp  |  -  3.44(  ^  fp\  =  ap[  -  3.44(  f\  (10) 
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representing  Fried's  (1966)  version  of  the  long  term  atmospheric  modulation  transfer 
function.  Equivalently,  Lutomirski  and  Yura's  (1971)  mutual  coherence  function  (MCF) 
is 

MCf(/>)-exp|-(-£-)5,J),  (11) 

where  p0  is  an  autocorrelation  length  of  the  electric  field. 

The  parameters  r0  and  p„  arc  different  measures  of  the  electromagnetic  (EM) 
wave's  coherence  as  restricted  by  the  atmospheric  turbulence.  They  are  related  by 

r0=  3.443/5p0  =  2.1p0.  (12) 


Using  Eq.  (8)  and  (9),  the  coherence  length  for  a  plane  wave,  r0  becomes 


n 


/•0  =  2.1[1.46A3 


cl(2)dzT,,s. 


•'O 


(13) 


The  parameter  r0  represents  the  distance  traverse  to  the  direction  of  propagation  where, 
on  average,  the  correlation  of  the  electric  field  is  e~yM  (Fried,  1966). 

To  interpret  r0  consider  the  following: 

1.  For  an  optical  aperture  d  smaller  than  the  coherence  length,  the  electric  field  across 
the  aperture  is  coherent.  The  angular  resolution  will  be  proportional  to  ).\d . 

2.  For  an  optical  aperture  larger  than  the  coherence  length,  the  correlation  across  the 
aperture  is  limited  to  regions  of  r0  or  p0  (depending  on  how  the  correlation  is  de¬ 
fined).  The  E.Y1  intensities  across  the  aperture  will  add  (energy  is  conserved).  The 
optical  system's  average  angular  resolution  will  be  proportional  to  /./r0  (Walters, 
Favier  and  Hines,  1979). 

3.  Visually,  for  apertures  smaller  than  r0,  atmospheric  turbulence  results  in  image 
centroid  motion,  a  process  called,  "tilt".  In  telescopes  larger  than  r„  the  image  is 
broken  up  into  multiple  speckels,  each  associated  with  a  r0  sized  region.  The  angle 
//r0  determines  the  overall  envelope  of  the  complex  image. 

4.  With  respect  to  measuring  the  atmospheric  optical  turbulence,  large  r0  values,  such 
as  200-400  mm,  indicate  a  small  amount  of  optical  turbulence  is  present  along  the 
integrated  path.  Conversely,  small  r0  magnitudes,  such  as  20-40  mm,  imply  a  large 
amount  of  optical  turbulence  exists  along  the  integrated  path. 

2.  Isoplanatic  Angle 

The  isoplanatic  angle,  00  ,  is  an  angular  measurement  of  spatial  coherence  be¬ 
tween  two  intersecting  rays  (Walters,  1985).  Given  two  light  sources  a  small  angle  0 
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apart,  the  random  atmospheric  phase  fluctuations  along  the  two  paths  are  different. 
Figure  2  displays  rays  r3,  and  r31  interacting  with  three  different  turbulence  patches.  At 
the  atmosphere  top  the  two  paths  undergo  unique  path  distortions.  At  the  surface  (near 
point  1),  the  turbulent  interaction  is  common  to  both  rays.  0O  is  the  angle  between  rays 
r2l  and  r31  such  that  the  accumulated  phase  fluctuations  are  correlated  to  within  e~'  of  a 
perfect  correlation. 

Fried  (1982)  develops  a  mathematical  representation  for  d0  for  apertures  d  >  > 
h  ■ 


0o  =  [2.91/c2 


Cl(2)2>l3dzfllS, 


(14) 


where  k  is  27://,  wavenumber  of  light;  Q  is  the  refractive  index  structure  parameter;  and, 
z  is  the  altitude  from  the  ground. 

Typical  Q  profiles  display  a  rapid  decrease  with  altitude.  The  z5'3  term  is  sig¬ 
nificant  in  that  it  emphasizes  the  optical  turbulence  at  higher  altitudes,  The  net  product 
of  these  two  terms  (the  integrand)  results  in  a  profile  of  relatively  constant  values  from 
a  few  kilometers  to  10  —  15  kilometers  (Walters,  1985). 

The  practical  interpretation  of  0O  is  similar  to  r0.  When  0O  magnitudes  arc  small 
(1-2  ju  rad),  the  optical  turbulence  along  the  integrated  path  and  especially  around  the 
tropopause  level  is  large.  Conversely,  when  0O  values  are  large  (14-20  n  rad),  optical 
turbulence  is  small. 


13 


Fig.  2.  Diagram  of  Isoplanalic  Angle,  0  0 . 
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IV.  EXPERIMENT/DATA 


%  A.  RME  KIHEI  DATA  COLLECTION 

The  initial  sampling  of  isoplanatic  angles,  60,  took  place  on  the  summit  of  Mt. 

*  Haleakala  in  Maui,  Hawaii.  These  measurements  commenced  on  17  September  1985 
and  concluded  on  1  February  1986.  Simultaneous  meteorological  surface  measurements 
at  the  site  provided  some  understanding  of  localized  atmospheric  events.  A  year  and  a 
half  later,  both  00  and  r0  were  measured  at  the  base  of  Haleakala  (at  approximately  sea 
level),  in  Kihei.  Six  sampling  sessions  of  four  and  five  nights  took  place  between  August 
19S7  and  July  19SS.  Two  Hawaii  Rawinsonde  Observation  stations,  Hilo  and  Lihue, 
provided  raw  data  listings  for  each  of  the  1987-1988  cases.  Additional  on-site 
rawinsonde  launches  supplemented  the  optical  data  for  four  out  of  the  six  sessions  (Ta¬ 
ble  1  ).  The  meteorological  data  for  the  August  and  December  1987  experiments  used 
the  VIZ  WO-SOOORP+  rawinsonde  system  with  10-sccond  vertical  resolution.  The 
April  and  July  1988  experiments  employed  the  Vaisala  DigiCORA  MW  11  system  with 

*■  2-  and  10-sccond  vertical  resolution.  All  four  cases  used  the  OMEGA  Navigational 

Network  for  wind  measurements.  Intermittent  GOES- WEST  Infrared  Satellite  Images 

♦  complimented  the  rawinsonde  data  for  all  the  cases.  The  primary  focus  of  this  study 
will  be  on  the  April  and  July  1988  experiments,  due  to  the  higher  vertical  resolution 
provided  by  the  DigiCORA  dataset. 


Table  1.  MAIN  OPTICAL  AND  METEOROLOGICAL  RME  SENSORS 
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B.  KIHEI  RME  EQUIPMENT 

1.  Optical  Equipment 

Optical  measurements  were  taken  by  an  isoplanometer  and  transverse  coherence 
length  sensor  designed/built  by  Dr.  D.L.  Walters,  NPS.  For  specific  details  on  each  in¬ 
strument  the  reader  should  pursue  Stevens  (19S5)  for  the  isoplanometer  and  Walters, 
Favier,  and  Hines  (1979)  for  the  transverse  coherence  length  sensor.  Isoplanatic  angles 
were  measured  by  an  isoplanometer  mounted  on  a  portable  200  mm  Celestron  telescope. 
The  software  was  run  on  an  HP  217  computer.  Isoplanometer  data  were  automatically 
sampled  once  every  second  and  displayed  real-time  on  a  time  verses  isoplanatic  angle 
plot.  Before  storing  the  data,  the  samples  were  converted  into  10  second  averages.  The 
complete  isoplanatic  angle  system  is  capable  of  operating  independent  of  human  inter¬ 
action  for  1-3  hours  provided  polar  alignment  is  perfect  and  the  chosen  star's  zenith 
angle  does  not  exceed  40  degrees. 

The  transverse  coherence  length  sensor  was  mounted  onto  a  portable  355  mm 
Celestron  telescope.  Software  used  for  real-time  processing  required  an  operator  to 
evaluate  the  image  and  stellar  positioning  before  each  sample.  This  added  interface  with 
a  trained  researcher  served  to  suppress  instrumental  artifacts.  The  average  sampling  rate 
was  between  1-2  minutes/sample.  Processed  data  were  stored  on  an  HP  300  series 
computer. 

2.  Meteorological  Equipment 

Meteorological  measurements  specific  to  this  thesis  were  a  product  of  the 
Vaisala  DigiCORA  MW  11.  This  solid  state,  digital  design  rawinsonde  package  inte¬ 
grates  the  receiver,  data  processors  and  operator  display  into  one  unit.  High  speed 
precision  sensors  permit  rapid  sampling  (at  approximately  2  second  intervals)  resulting 
in  high  density  data  acquisition.  The  95  mm  by  145  mm  (approximate  size)  RS  80 
Rawinsonde  (Navaid  Radiosonde)  package  was  launched  on  an  unravelling  tether  at¬ 
tached  to  a  200  gm  balloon  filled  with  helium.  Pressure  values  were  measured  by  a 
BAROCAP  pressure  sensor,  a  steel  capacitive  aneroid  with  a  resolution  of  0.1  hPa  and 
an  accuracy  of  0.5  hPa.  Temperature  sampling  was  taken  by  a  THERMOCAP  sensor, 
a  small  capacitive  bead  in  a  glass  casing.  Resolution  of  the  THERMOCAP  tvas  0.1C; 
accuracy  was  0.2C.  The  H  CM  I  CAP,  a  thin  film  humidity  sensor,  measured  humidity 
with  a  resolution  of  1%  relative  humidity  (RI  I)  and  an  accuracy  of  2%  RH.  According 
to  the  Vaisala  Ref.  A0435,  this  sensor  produced  a  reliable  response  even  in  low  temper¬ 
atures  and  after  exposure  to  condensation.  The  measuring  range  published  was  0-100% 
RH  (Vaisala,  1986a). 
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The  commutation  between  the  sensors  and  the  transducer  electronics  was  ac¬ 
complished  by  a  digital  electronic  switch.  A  water  activated  battery  provided  the  RS  80 
"sonde"  with  a  light  weight  power  supply  capable  of  keeping  the  package  activated  for 
the  entire  surface  to  30  km  (balloon  permitting)  launch  sequence  (Vaisala,  1986b). 

The  telemetry  frequency  was  between  400-406  MHz.  Wind  data  for  this  exper¬ 
iment  were  calculated  from  the  OMEGA  Navigational  System.  The  RS  80  received 
signals  from  the  OMEGA  Navaid  Network  (11-14  kHz),  superimposed  this  information 
onto  the  thermodynamic  data  and  relayed  the  sampling  to  the  DigiCORA  MW  11 
(Vaisala,  1986b). 

C.  KIHEI  SONDE  DATA  PREPARATION  FOR  ANALYSIS 

The  initial  form  of  the  DigiCORA  MW  11  rawinsonde  output  used  in  this  research 
began  as  a  hardcopy  listing  of  time  in  minutes  and  seconds,  height  in  meters,  pressure 
in  millibars,  temperature  in  Celsius,  relative  humidity  in  percent,  dewpoint  temperature 
in  Celsius,  mixing  ratio  in  grams  per  kilogram,  wind  speed  in  knots  and  wind  direction 
in  degrees.  Using  a  Dest  PCSCAN  scanner,  the  powerful  IBM  mainframe  XED1T,  a 
software  package  capable  of  transfering  IBM  files  to  a  IBM  PC- AT,  and  the  instant 
graphics  display  of  QUATTRO,  the  initial  hardcopy  was  duplicated  as  an  ASCII  file. 

Thermodynamic  sampling  rate  for  the  first  two  kilometers  was  every  two  seconds. 
Beyond  this  threshold,  sampling  continued  once  every  10  seconds.  Due  to  the  limita¬ 
tions  in  the  Omega  Navigational  system,  winds  were  measured  every  10  seconds 
throughout  the  launch.  A  linear  interpolation  between  the  10-second  wind  measure¬ 
ments  was  applied  to  the  initial  two  kilometers  remedying  the  inconsistent  sampling  rate. 

Variables  calculated  from  the  raw  data  include  potential  temperature,  potential 
temperature  gradient,  shear,  Richardson  number,  Q- ,  Q,  >•„  and  0o . 
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V.  ANALYSIS/RESULTS 


To  strengthen  the  link  between  the  measured  optical  parameters  and  the 
meteorological  data  collected,  a  number  of  methods  for  calculating  r0  and  60  values  were 
pursued.  As  explained  earlier,  both  r0  and  0O  are  a  function  of  Q  ,  which,  in  turn,  is  a 
function  of  Cf  (see  Chapter  3).  In  this  chapter,  four  methods  for  calculating  Cf  from 
rawinsonde  data  will  be  investigated.  However,  the  first  step  is  to  define  an  appropriate 
conservative  temperature  quantity  to  be  used  in  the  Cf  calculation. 

A.  STANDARDIZING  THE  TEMPERATURE  PROFILE 

An  air  parcel  vertically  displaced  in  the  atmosphere  experiences  a  temperature 
change  as  the  parcel's  pressure  equalizes  with  the  ambient  environment  (Clifford,  1978). 
Permitting  the  air  to  be  moist  but  unsaturated,  and  assuming  the  atmosphere  is  in 
hydrostatic  equilibrium,  the  rate  at  which  the  lifted  parcel's  temperature  will  decrease  is 
the  dry-adiabatic  lapse  rate.  This  is  approximately  9.767C.’km.  If  the  parcel  is  saturated, 
the  release  of  latent  heat  produces  a  lapse  rate  that  is  less  than  the  dry-adiabatic  lapse 
rate  (Huschkc,  1970).  For  this  study,  it  is  important  to  work  with  conservative  variables, 
a  variable  that  is  independent  of  position  in  space.  As  described  above,  temperature 
does  not  fit  this  description.  The  conservative  quantity  (independent  of  position)  best 
suited  is  potential  temperature  (Clifford,  1978).  By  definition,  potential  temperature  is 
"the  temperature  which  the  parcel  of  air  would  have  if  it  were  expanded  or  compressed 
adiabatically  from  its  existing  pressure  and  temperature  to  a  standard  pressure  p0  (gen¬ 
erally  taken  as  lOOOmb)"  (Wallace  and  Hobbs,  1977).  In  short,  the  temperatures  sam¬ 
pled  vertically  are  standardized  to  one  pressure  level.  Mathematically,  the  potential 
temperature  expression  is 


(15) 

where  0  is  potential  temperature  in  kclvin,  T  is  temperature  in  kelvin,  P0  is  a  standard 
pressure  such  as  lOOOmb,  P  is  pressure  in  millibars,  R  is  the  gas  constant  for  dry  air  and 
Cf  is  the  specific  heat  of  dry  air  at  constant  pressure  (1.005  J/g  K)  (Husehke,  1970). 
Note:  Following  the  usual  convention,  Q  shall  hereafter  be  written  in  place  of  C| 
(Clifford,  1978). 
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B.  CALCULATING  Clr :  FOUR  CASE  STUDIES 

1.  Case  I:  Simple  Potential  Temperature  Differences 

Case  I  takes  the  given  mathematical  definition, 

,1  <1(0!  -§>)  — <®|— 5,)f> 

CT - ^3 -  <16) 

and  assumes  02  =  ©j  .  Physically,  this  presumes  homogeneity  between  neighboring 
samples.  In  other  words,  the  samples  are  taken  within  the  inner  and  outerscale  of  the 
turbulent  layer  (Tatarski,  1961).  Note:  the  angled  brackets  indicate  an  ensemble  average 
of  the  enclosed  quantity.  Associated  with  the  above  and  all  the  cases  hereafter  is  the 
assumption  that  the  sampled  layer  is  isotropic  (Tatarski,  1961). 

Unlike  most  Ot  calculations,  the  dataset  employed  by  this  research  consists  of 
vertical  temperature  samples.  For  the  above  assumption  to  be  true,  the  rate  of  sampling 
must  be  sufficiently  high  with  respect  to  the  potential  temperature  gradients.  If  the 
sampling  rate  is  too  slow,  trends  in  potential  temperature  will  show  up  in  the  differences, 
(02  -  0J.  Figure  3  displays  such  trends  first  between  2  and  5  km  then  again  above  12 
km.  Note:  The  large  increase  in  potential  temperature  differences  in  the  stratosphere 
was  found  in  all  the  profiles  examined. 

a.  Results 

Using  the  entire  potential  temperature  profile  for  the  Ot  calculations,  the 
eight  April  1988-K.ihei  rawinsonde  flights  and  the  ten  July  1988-Kihei  launches  all  under 
estimate  the  direct  /•„  and  60  measurements  by  a  factor  of  ten.  (See  Appendix  A)  This 
is  expected  considering  the  trends  in  the  potential  temperature  differences.  As  demon¬ 
strated  in  Figure  3,  these  trends  generate  differences  approaching  2  K  near  the  boundary 
layer  top.  In  the  stratosphere  (above  15  km),  potential  temperature  differences  are  even 
larger.  For  this  particular  profile,  the  stratospheric  potential  temperature  differences 
reach  as  much  as  4  K.  When  these  large  differences  are  blindly  used  to  calculate  Q-,  the 
result  is  a  deceptively  large  over  estimation  of  Ot  .  The  over  estimated  Q-  subsequently 
leads  to  an  under  estimation  of  r0  and  0a. 

2.  Case  II:  Using  a  Running  Mean 

Case  II  investigates  the  effects  of  detrending  the  potential  temperature  profile 
prior  to  calculating  Ot.  The  two  representative  high  pass  filters  selected  for  this  case 
include  the  equally  weighted  111  Running  Mean  (RM) 
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and  the  weighted  13531  Running  Mean 


A/=  .  U°J 

Figure  4  displays  the  potential  temperature  difference  profile  after  the  1 1 1  RM  has  been 
subtracted  from  the  original  potential  temperature  profile.  There  is  still  an  increase  in 
potential  temperature  differences  around  the  boundary  layer,  however,  the  peak  values 
of  approximately  0.3  K  is  balanced  with  an  adjacent  value  of  about  -0.3  K.  Even  the 
stratosphere  appears  nearly  symmetrical  between  +/•  0.7  K. 

a.  Ill  Running  Mean 

The  implementation  of  the  equally  weighted  (111)  RM  handles  the 
DigiCORA's  two-  to  ten-second  sampling  rate  transition  by  using  values  equally  spaced 
in  time  when  available.  When  data  are  unavailable,  preference  is  given  to  the  "central" 
point.  Profile  endpoints  arc  ignored.  For  two-second  data  the  filter  takes  into  account 
layers  between  15-30  meters  thick.  Ten-second  data  involves  layers  of  90-120  meters. 

b.  13531  Running  Mean 

The  13531  running  mean  filter  accommodates  the  two-  to  ten-second  sam¬ 
pling  rate  change  with  preference  given  first  to  data  equally  spaced  in  time,  and  second, 
to  the  central  point.  Endpoints  are  noted,  but  generally  ignored. 

c.  r0  Results 

Both  RM  methods  produce  r0  values  within  the  same  order  magnitude  as 
the  optical  sensor  values.  With  few  exceptions,  the  relatively  large  and  small  magnitude 
variations  present  in  the  optical  sensor  data  are  traced  by  both  RM  results.  The  13531 
RM  generate  values  approximately  16%  larger  than  the  1 1 1  RM.  The  error  between  the 
111  RM  and  sensor  results  ranges  from  13%  to  55%.  The  13531  RM  under  estimates 
the  sensor  values  by  4%  to  48%.  When  the  Q-  profile  detrended  by  the  111  RM  is 
truncated  just  below  the  stratosphere's  large  Ot  increase,  r0  increases  by  Only  11%.  Ta¬ 
ble  2  displays  a  sample  of  these  results. 

d.  0  o  Results 

While  large  scale  magnitude  changes  recorded  by  the  sensor  are  virtually 
ignored,  the  calculated  0j  values  are  within  the  same  order  of  magnitude  as  the  sensor. 
The  range  of  sensor  measurements  is  from  3  to  7.5  n  rad.  Values  calculated  from  the 


Ill  RM  cover  a  range  of  0.7-1.3  n  rad.  (The  8.5  n  rad  of  flight  5,  July  88  is  ignored.) 
0O  calculated  with  the  13531  RM  ranges  from  0.8-1.5  n  rad.  Truncating  the  detrended 
Ot  profile  just  below  the  stratosphere's  large  Ot  increase,  increases  80  by  78%.  A  sample 
of  these  results  is  in  Table  2  . 


Table  2.  SINGLE  PROBE  RESULTS  -  0605  UTC  28  APRIL  1988 


Or  CALCULATION 
METHOD* 

r0  (mm) 

0„  (M  rad) 

SPTD  (all  values) 

8 

0.2 

BG  (all  values) 

7 

0.2 

1 1 1  RM  (all  values) 

36 

1.3 

111  RM 

(troposphere  only) 

40 

2.3 

13531  RM  (all  values) 

3S 

1.5 

SENSOR  VALUES 

43 

3 

3.  Case  III:  Using  the  Bulk  Gradient  C>  Equation 

The  third  approach  incorporates  an  expression  developed  by  Tatarski  (1961) 
that  relates  Or  to  bulk  parameters  of  the  atmosphere.  This  expression  which  will  here¬ 
after  be  referenced  as  the  "Bulk  Gradient  Cf "  equation  (BG)  is, 

c\=<mLAJ\dQldzf,  (19) 

where  a  is  a  constant;  a  is  the  ratio  of  eddy  thermal  difiusivity  to  eddy  viscosity;  L0  is  the 
outerscale  length;  and.  (dQjdzy  is  the  potential  temperature  gradient.  Because  L0  is  not 
universally  well-defined  and  dQjdz  is  also  vague  with  respect  to  the  resolution  needed  to 
properly  define  the  derivative,  the  values  calculated  from  BG  are  suspect.  The  main 
strength  of  the  expression  is  that  it  accentuates  the  strong  dependence  Cf  has  on  the 
square  of  the  temperature  gradient.  This  term  also  provides  the  "atmospheric  density 
inhomogeneity  susceptible  to  being  mixed  by  the  wind  turbulence"  (Walters  and  Kunkel, 
1981).  Assumptions  made  when  applying  this  equation  include  the  following: 

1.  The  atmosphere  is  mechanically  turbulent.  Note  the  lack  of  a  shear  term. 

2.  The  outerscale,  7„0,  is  half  the  distance  between  samples. 

3.  Constant,  a,  is  2.S. 

4.  Ratio,  a  ,  is  1  . 
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a.  Results 

The  BG  r0  and  0O  results  are  one-tenth  the  sensor  measurements.  Both  op¬ 
tical  parameters  calculated  in  this  case  closely  trace  those  values  generated  by  the  Simple 
Potential  Temperature  Differences  (SPTD)  Case.  r0  magnitudes  vary  over  a  range  of  6-8 
mm.  0O  values  cover  a  range  of  2-3  mm.  These  values  are  within  7%  of  the  SPTD  case. 
Considering  the  dominant  term  of  the  BG  equation  is  the  potential  temperature  gradient, 
it  is  not  surprising  to  discover  the  final  results  are  so  closely  associated  to  Case  I.  An 
explanation  for  why  r0  and’0o  values  are  off  by  a  factor  of  10  begins  with  a  re-evaluation 
of  the  assumptions.  The  constant,  a,  has  a  suggested  range  between  2.8  to  3.2 .  Shifting 
this  to  the  higher  constant  will  not  accommodate  the  discrepancy.  Changing  the 
outerscale  length  to  the  full  thickness  of  the  layer  serves  only  to  decrease  r0  and  00  . 
Justifying  an  outerscale  smaller  than  half  the  thickness  of  the  layer  sampled  taps  into  the 
fundamental  problem  with  the  Ot  calculations.  That  is,  vertical  resolution.  The  sampled 
values  are  too  far  apart  to  provide  the  temperature  gradient  information  needed  to  drive 
this  CJr  equation  properly.  Decreasing  the  outerscale  length  without  simultaneously  de¬ 
creasing  the  thickness  of  the  sampled  layer  is  illogical.  As  in  Case  1,  the  BG  results  in¬ 
dicate  a  need  for  either  detrending  the  data  prior  to  calculating  Ot  or  using  a  rawinsonde 
dataset  capable  of  producing  more  detailed  temperature  gradient  information. 

4.  Case  IV:  Using  Richardson  Number 

Cf  is  a  constant  representing  the  total  amount  of  energy  in  a  turbulent  layer 
(Clifford,  1978).  The  vertical  temperature  samples  measured  by  the  rawinsonde  are  not 
always  taken  from  turbulent  layers.  To  computationally  isolate  only  those  atmospheric 
layers  which  are  turbulent,  the  Richardson  number,  Ri,  wfas  used.  Richardson  number 
is  "the  ratio  of  work  done  against  gravitational  stability  to  energy  transferred  from  mean 
to  turbulent  motion"  (Huschke,  1970).  When  Ri  <  0.25,  the  layer  is  unstable  and  Cf  is 
calculated  using  the  SPTD,  Running  Mean  and  BG  approaches,  respectively.  When  Ri 
2,  the  layer  is  stable  and  no  C2r  calculation  is  made.  "Critical"  conditions  exist  when 
0.25  <  Ri  <  2.0  .  As  the  final  case  reviewed,  Ot  is  calculated  when  Ri  <  2.0  . 
a.  Results:  Unstable  Ri  Case  (Ri  <  0.25) 

Letting  a  turbulent  layer  be  defined  as  one  having  a  Ri  <  0.25,  each  of  the 
above  three  methods  for  calculating  Ot  w'as  conditionally  computed.  With  remarkable 
consistency,  all  three  methods  over  estimate  r0 .  r0  values  farthest  from  the  sensor's  value 
generally  are  produced  by  the  SPTD  case.  Except  for  one  extreme  case  where  r0  is  cal¬ 
culated  as  greater  than  360  mm  from  the  sensor  value,  the  SPTD  r0  values  are  generally 
2  to  3  times  too  large.  The  BG  method  leads  to  values  that  are  approximately  7%  closer 
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than  SPTD  to  the  sensor  value.  The  RM  method(s)  r0  values  over  estimate  the  sensor 
by  18  to  71  mm  for  the  13531  RM  and  8  to  58  mm  for  the  111  RM. 

0O  calculations  generate  less  consistent  results.  The  sequence  of  greater  to 
smaller  magnitudes  seen  in  the  r0  results  does  not  hold  for  0o .  About  40%  of  all  the  0o 
data  calculated  for  Case  IV  (Ri<0.25)  over  estimate  the  sensor's  measurement;  50% 
under  estimate  0O;  and,  10%  straddle  the  sen^r  value  (see  Appendix  A) . 
b.  Results:  Conditionally  Unstable  Ri  Case  ( Ri  <  2.0) 

With  the  inclusion  of  more  atmospheric  layers  in  the  r0  integral,  one  would 
expect  the  resulting  magnitude  to  be  smaller.  For  the  most  part  this  is  true.  The  SPTD 
and  the  BG  cases  respond  in  the  same  manner.  Their  r0  values  under  estimate  the  sensor 
by  17-80%  for  the  SPTD  case  and  22-65%  for  the  BG  case.  The  111  and  13531  RM 
calculations  still  over  estimate  r0  with  the  13531  RM  slightly  higher  than  the  111  RM 
(See  Appendix  A).  The  magnitude  of  the  results  as  compared  to  the  sensor's  value  is 
generally  within  78%  for  13531  RM  and  59%  for  the  111  RM. 

0 o  results  display  a  pattern  consistent  with  the  "All  Values"  results.  The 
largest  to  smallest  00  magnitudes  calculated  are:  (Sensor  Value),  13531  RM,  111  RM, 
SPTD  and  BG.  All  values  are  less  than  3  n  rad!  As  before  the  large  magnitude  vari¬ 
ations  observed  in  the  sensor  data  arc  virtually  ignored  in  the  calculated  data. 

C.  CORRELATION  OF  SENSOR  AND  CALCULATED  PARAMETERS 

Correlating  the  measured  and  calculated  optical  parameters  for  each  case  discussed 
in  section- B  generates  r-squared  values  between  0.0009  and  0.5  .  Table  3  lists  r-squared 
values  for  each  case.  Note:  r-squared  is  1  for  a  perfect  correlation. 

The  lack  of  good  correlation  demands  a  review  of  the  sensor's  ability  to  sample 
values  for  Ot  calculations.  After  detrending  the  potential  temperature  profile  with  a  1 1 1 
RM,  potential  temperature  differences  for  the  two-second  data  are  plotted.  (See 
Figure  5)  Initially,  this  profile  appears  to  contain  valid  turbulence  information.  Re¬ 
peating  this  procedure  with  the  temperature  (K)  profile  (Figure  6),  however,  unveils  a 
major  problem.  The  systematic  +/-0.03  K  oscillation  dominating  the  profile  implies 
that  the  temperature  differences  below  2  km  are  predominantly  quantization  noise  from 
the  0.1°C  temperature  resolution.  Comparing  Figure  5  and  Figure  6sho\vs  that  the 
potential  temperature  calculation  introduces  additional  errors  that  obscure  the  quantiz¬ 
ing  error.  Rather  than  pursue  a  detailed  correlation  analysis,  attention  is  diverted  to  1) 
identifying  and  quantizing  the  major  sources  of  error,  and  2)  providing  recommendations 
for  a  temperature  sensor/system  more  suited  for  calculating  Q- . 
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Fig.  5.  Initial  2  km  of  Detrended  Potential  Temperature  Differences 


26 


Table  3.  MEASURED  AND  CALCULATED  PARAMETERS  CORRELATED 


R-SQUARED 

VALUES 

r0  (mm) 

6  o  (m  rad) 

C\  CALCU¬ 
LATION 
METHOD 

ALL 

Values 

Ri<0.25 

Ri  <  2.0 

ALL 

Values 

Ri<0.25 

Ri  <  2.0 

SPTD 

0.009 

0.3 

0.5 

0.05 

0.3 

0.02 

BG 

0.009 

0.3 

0.4 

0.05 

0.3 

0.02 

111  RM 

0.005 

0.0009 

0.06 

0.01 

0.2 

0.1 

13531  RM 

0.002 

0.001 

0.07 

0.07 

0.2 

0.2 

D.  ERROR  ANALYSIS 

1.  Temperature  Probe  Time  Constant  Correction 

In  both  the  April  SS  and  July  88  soundings,  the  temperature  was  measured  with 
a  capacitive  bead.  For  sensor  specifications  see  Chapter  4.  To  calculate  the  sensor  re¬ 
sponse  time  t  ,  the  following  equation  was  used: 

l/r  =  /l  +  Z?/v,  (20) 

where  v  is  velocity,  A  is  the  constant  for  the  calm  atmosphere  (v  =  0)  and  B  is  a  constant 
at  velocity,  v  .  For  this  particular  probe,  the  sensor  time  constant  at  a  velocity  of  6  m/s 
is  2.5  s.  With  an  ascent  rate  of  4  m/s,  t  increases  to  approximately  3  s.  In  order  to 
calculate  the  error  induced  by  the  temperature  sensor  time  constant,  the  structure  func¬ 
tion  of  the  probe  system  with  limiting  scale  lengths  must  be  compared  to  the  structure 
function  over  all  scale  lengths. 

Mathematically,  the  structure  function  over  all  scale  lengths  is  (Tatarski,  1961): 


D(y,0,oo)  =  4 


sin2{krl2)r5lictk, 


Jo 


(21) 


where  k  is  2  n  !). ,  /  the  turbulent  wavelength  and  r  the  probe  separation  distance. 

The  structure  function  filtered  by  a  one  pole,  low  pass  filter  is  (Huntcn,  1964): 


2S 


(22) 


rb 


D{r,a,b)  -  4 


*>a 


sin2(Ar/2)A:~5/3 
1  +  (wt)2 


dk, 


where 

a  is  2 it/Lo,  the  limiting  lower  frequency  (outerscale,  Lq), 

b  is  27t//0,  the  limiting  upper  frequency  (inner  scale,  /0), 

to  is  kV  where  V  is  the  wind  velocity  (rawinsonde's  ascent  rate)  and 

r  is  the  time  constant  of  the  temperature  probe  (the  upper  frequency  cutoff). 

The  cur  term  in  the  denominator  attenuates  the  high  spatial  frequencies.  Since  the  k~s'3 
term  increases  at  low  wavenumber,  Eq.  (22)  is  much  more  sensitive  to  the  selection  of 
an  outerscale  than  the  innerscale.  Test  calculations  of  the  D(8m,a,b)  integral  confirm 
that  an  innerscale  of  0.1  m  verses  0.001  m  produce  virtually  identical  results.  When 
varying  the  outerscale  from  10  m  to  1000  m  (all  other  parameters  of  D(8m,a,b)  remain¬ 
ing  unchanged),  the  1000  m  structure  function  is  200  times  larger  than  for  the  10  m. 

With  an  inner  scale  of  0.1  m,  an  outerscale  of  1000  m  and  a  3-second  temper¬ 
ature  probe  time  constant,  the  DigiCORA's  2-second  data  is  15%  of  the  actual  Ot 
present.  The  structure  parameter  Q,  calculated  from  the  10-second  data  is  54%  of  the 
actual  Ot  value.  While  the  10-second  data  appears  to  produce  results  more  than  three 
times  closer  to  the  desired  values,  this  is  somewhat  deceptive.  The  larger  time  interval 
between  samples  (accent  rate  is  constant)  results  in  a  greater  probe  separation  distance. 
Consequently,  higher  spatial  frequencies  or  rapid  changes  in  the  potential  temperature 
will  be  missed.  This  aliasing  of  temperature  measurements  is  one  of  the  major  difficulties 
with  calculating  structure  function  parameters  from  rawinsonde  data. 

Another  problem  with  calculating  a  structure  function  D  from  10  second  (r=40 
m)  data  involves  a  limitation  of  the  Kolmogorov  rm  slope  model.  D  is  a  function  of  the 
probe  separation  distance,  r,  and  the  outerscale.  For  r  values  comparable  to  the 
outerscale,  the  actual  D(r)  deviates  from  Kolmogorov's  model.  When  r  is  larger  than 
the  outerscale,  the  measured  (actual)  structure  function  will  be  smaller  than  the  D  cal¬ 
culated  from  Kolmogorov's  model.  (See  Figure  7)  The  significance  of  this  problem  is 
that  the  calculated  D  under  estimates  the  actual  turbulence  present  when  r  is  greater 
than  the  outerscale.  Thus,  accepting  the  calculated  Cf  as  54%  of  the  actual  Cf  present 
is  done  cautiously. 
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a.  Application  of  the  Filter 

Correcting  the  basic  Q  equations  for  the  temperature  probe  time  constant 
produces  an  even  greater  over  estimation  of  atmospheric  turbulence.  The  magnitudes 
of  the  filtered  r0  and  0O  (all  values)  decrease  35-50%  from  the  unfiltered  values. 
Table  4  and  Table  5  display  the  results  taken  from  two  rawinsonde  datasets:  0605  UTC 
28  April  1988  and  1434  UTC  30  April  1988.  Imposing  this  correction  onto  the  calcu¬ 
lation  of  the  optical  parameters  where  Ri  <  0.25  and  Ri  <  2.0  produce  the  same  35-50% 
decrease  as  the  "All  Values"  cases.  Unless  otherwise  noted,  further  reference  of  cor¬ 
rection  implementation  will  be  done  on  the  0605  UTC  28  April  1988  Kihei,  Maui 
rawinsonde  dataset. 


Table  4.  CORRECTING  FOR  TEMPERATURE  PROBE  TIME  CONSTANT 


0605  UTC  28  April  1988 

r0  (mm) 

0 o (m  rad) 

Q-  Calculation  Method 

Corrected 

Uncor¬ 

rected 

Corrected 

SPTD-A11  Values 

8 

5 

0.2 

0.1 

BG-A11  Values 

7 

5 

0.2 

0.1 

111  RM -All  Values 

36 

20 

1.3 

0.9 

13531  RM-A11  Values 

38 

22 

1.5 

1.0 

SENSOR  VALUES 

48 

48 

3 

•y 

Table  5.  CORRECTING  FOR  TEMPERATURE  PROBE  TIME  CONSTANT 


1434  UTC  30  April  1988 

>o  (mm) 

0O  (m  rad) 

Ot  Calculation  Method 

Uncor¬ 

rected 

Corrected 

Uncor¬ 

rected 

Corrected 

SPTD-All  Values 

8 

4 

0.2 

0.1 

BG-A11  Values 

8 

4 

0.2 

0.1 

111  RM-A11  Values 

33 

16 

1.2 

0.8 

13531  RM-A11  Values 

37 

19 

1.3 

0.9 

OPTICAL  SENSOR 
VALUES 

54 

54 

7 

7 

31 


2.  Sensor/ Atmosphere  Interaction  Noise 

Since  the  optical  parameters  calculated  by  the  rawinsonde  are  smaller  than  the 
optical  sensor  values  and  compensation  for  the  probe  response-time  accentuates  this 
problem,  some  other  large  error  source  must  be  present.  Noise  introduced  by  the 
sensor-atmosphere  interaction  is  such  a  source.  Two  methods  attempt  to  quantize  this 
phenomena. 

a .  Standard  Error  Formula 

The  first  approach  tackles  the  problem  by  focusing  on  the  limitations  of  the 
rawinsonde  system.  In  other  words,  knowing  the  resolution  of  each  variable  measured, 
how  well  can  the  given  instruments  measure  atmospheric  turbulence?  This  method  uses 
the  standard  error  formula  which  for  a  function  E(x,y,z)  is  (Bevington,  1969): 


4=( 


0 E  \2/  \2  |  / 

ar)W  +( 


■f-)V  +  (-§7)W 


(23) 


Applying  the  standard  error  formula  to  the  SPTD  calculation  of  Ot  produces 


oil  =  4 
cr 


(T,  -  r,) 
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Note:  T2  and  7,  are  assumed  to  generate  independent  random  errors. 
With  Ot  defined  as 


Cn  —  (79  x  10-6  fcj-. 


the  system  induced  error  for  Q  is 
(79  x  icr6)V 


2 

acr 


(PVy+(if^r^+(2c^ 


(25) 


(26) 


Implementing  the  error  correction  reveals  three  layers  where  (Or-  ocfi  <  0.  These  re¬ 
sults  imply  that  the  measurement  uncertainties  are  comparable  to  the  measurements. 
The  sensitive  regions  include  0  to  2  km,  4.5  to  9.5  km  and  12  to  14  km.  However,  the 
area  of  major  concern  is  near  the  surface  where  the  majority  of  turbulence  is  measured. 
Eliminating  the  0-2  km  layer  would  result  in  a  gross  under  estimation  of  turbulence 
which  ultimately  generates  large  r0  and  0O .  Since  by  definition  Cf  can  not  be  less  than 
zero,  all  negative  error  corrections  are  assumed  to  be  zero. 
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The  results  from  the  Q  error  calculation  show  values  less  than  or  equal  to 
0  near  the  surface  (0  -  2  km).  Between  3.5  -  9.5  km,  as  well  as,  10  -  14  km,  the  Q  <  0 
condition  is  intermittent.  Unlike  the  Ot  calculation,  Cjj  has  a  pressure  term  that  effec¬ 
tively  minimizes  the  influence  of  any  values  above  the  tropopause.  Thus,  the  main 
concern  is  with  turbulence  measurements  near  the  surface. 

The  average  Ot  structure  parameter,  including  sensor  error  (0605  UTC  28 
April  1988),  is  (1.1  +  1.8)jcl0_3/?wr2/3  for  the  0  -  2  km  two-second  data;  the  average 
Q  ±  sensor  error  over  this  same  region  is  (8  ±  13)*10-16m-2''3 .  The  10-second  data  (2  km 
and  above)  produces  an  average  Ot±  sensoi  error  of  (7.4  ±  0.8)*10'2A?m"2'3;  the  average 
Q,  including  sensor  error  for  the  10-second  data,  is  (2.0  ±  0.4).rl0-’5m-2/3.  Clearly,  the 
0-2  km  data  are  buried  in  uncertainties. 

b.  Vertical  C>  Quantization  Noise  Simulation 

The  second  technique  used  for  quantizing  the  roundoff  error  in  the  SPTD 
(Simple  Potential  Temperature  Differences)  calculation  of  Ot  simulates  the  atmosphere's 
temperature  profile  with  a  straight  line  slope  of  -0.007  K/m.  This  technique  (Walters, 
1989)  includes  a  correlation  between  adjacent  temperature  samples  that  is  ignored  in  the 
standard  error  approach.  The  0.1  K  temperature  resolution  produces  a  quantization 
noise  of  1.4.\i0'3A'2/7r2/3  for  the  two-second  data  and  6.S.vl0~3A’2/7i_2?3  for  the  ten-second 
data.  The  two-second  noise  is  22%  less  than  the  Standard  Error  Function's  result.  The 
ten-second  value  is  the  same.  Running  the  simulation  on  a  detrended  potential  temper¬ 
ature  profile,  the  calculated  noise  for  Q-  is  a  factor  of  10  smaller  than  the  SPTD  calcu¬ 
lation.  Specifically,  the  two-  and  ten-second  quantization  noises  are  6.2*1 0_4A*w2/3  and 
.  1.4*10'4A2»r2'3,  respectively. 

3.  Implementing  Noise  and  Filter  Corrections 
The  measured  Q  is 


c). 

•  r 


mtosurtd 


(Filter  x  c\  ,)  +  Noise, 

x  *  actual'  1 


(27) 


where  the  Filter  is  the  temperature  probe  time  constant  correction,  and  Noise  is  the 
vertical  Ot  quantization  noise  from  the  simulation. 

Applying  the  noise  correction  to  the  111  RM  Ot  calculation  for  the  0605  UTC 
28  April  1988  rawinsonde  flight,  the  original  36  mm  r0  calculated  by  this  method  in¬ 
creases  to  46  mm  (all  values  included).  The  sensor  value  is  48  mm.  Replacing  all  neg¬ 
ative  Q-  with  zero  reduces  this  measurement  to  41  mm.  Imposing  the  temperature  probe 
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time  constant  correction  reduces  this  value  again  to  36  mm  (all  values  included)  and  ul¬ 
timately  r„  is  25  mm  including  only  positive  C?  values. 

60  displays  a  much  smaller  range  of  results  due  to  the  equation's  sensitivity  for 
values  at  the  tropopause  coupled  with  a  lack  of  large  variation  created  by  the  noise  and 
filter  at  that  level.  Specifically,  0O  starts  as  1.3  n  rad  (only  the  111  RM  applied).  In¬ 
cluding  all  values  and  subtracting  out  the  quantization  noise  increases  0O  to  1.4  n  rad. 
Removing  all  negative  Cf  does  not  change  the  calculated  measurement.  Dividing  the 
noise-free  Cf  by  the  filter  leads  to  a  value  of  1.0  n  rad  (all  values).  And  finally,  using 
only  positive  Cf,  the  60  result  remains  1.0  n  rad.  The  sensor's  measurement  is  3  n  rad. 

E.  DETERMINING  AN  ADEQUATE  TEMPERATURE  RESOLUTION 

The  (0. 1  ±  0.05 )K  temperature  probe  resolution  dominates  the  error  term  in  Eq. 
(24).  An  important  consideration  is  to  determine  the  temperature  resolution  required  for 
calculating  Cf  from  rawinsonde  data.  First,  let  <rcj.=  I.vIO^Attt2'3 .  A  Cf  noise  level  of 
l.YlO-^Am-2''3  is  adequate  for  most  measurements.  Assuming  the  temperature  difference 
is  0.0098  K  over  a  I  meter  height  difference  (dry-adiabatic  lapse  rate),  and  that 
on  =  =  0,  one  can  solve  for  of.  Note:  of2  =  a\x.  The  result  is  c\—  1.30.Yl0'5A?«r2'3. 

The  roundoff  error  (noise)  temperature  variance  expression  (Bose,  1985) 

4— iT!0"M.  (») 

where  d  is  the  number  of  base  ten  significant  digits  to  the  right  of  the  decimal  point  re¬ 
maining  after  roundoff,  permits  the  calculation  of  significant  digits  required  for  a  given 
variance.  Solving  for  d,  the  resolution  necessary  for  a  ocj.=  LvlO‘4Awr2'3  is  approxi¬ 
mately  two  decimal  points  remaining  after  roundoff,  0.01  K  resolution.  Repeating  the 
above  procedure  with  a  ATjAz  of  -0.007K,'lm  (a  lapse  rate  observed  in  the  Kihci  data¬ 
set),  the  results  are  also  approximately  two. 

The  initial  use  of  the  standard  error  formula  defined  aT  as  0.05  K,  as  implied  by  the 
sensor  specifications.  The  roundoff  error  temperature  variance  expression  calculates  oT 
as  0.03  K.  When  0.03  K  is  used  in  the  Cf  standard  error  formula,  the  product  is  ap¬ 
proximately  one-third  the  magnitude  of  the  standard  error  formula  using  0.05  K.  This 
decrease  in  the  Cf  noise  term  increases  Cf.  A  larger  Cf  leads  to  a  larger  Q  and,  subse¬ 
quently  a  smaller  /•„  and  0O.  Since  the  calculated  optical  parameters  are  already  under 
estimating  the  sensor  values,  this  adjustment  only  confirms  that  the  current  temperature 
sensor/system  is  at  best  marginally  acceptable  for  calculating  Cf  . 
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F.  RUNNING  MEAN  FITTING  ERROR 

A  fitting  error  around  sharp  gradients,  such  as  across  inversions,  arises  when  im¬ 
plementing  the  running  mean  into  the  C5r  calculation.  The  running  mean  is  unable  to 
follow  the  trends.  Figure  8  displays  the  large  difference  that  occurs  between  the  running 
mean  and  a  given  profile.  Though  variances  are  used  in  the  (^calculation,  the  variances 
are  artificially  enlarged  by  the  poor  fit  across  the  inversion.  The  SPTD  experiences  a 
form  of  this  fitting  problem  that  results  in  aliasing.  (See  Section  D1  of  this  chapter). 


Pressure— Altitude  (km) 


0605  UTC  28  April  80 


Fig.  8.  1 11  Running  Mean  Across  an  Inversion 


VI.  DISCUSSION  AND  CONCLUSIONS 


A.  CORRELATED  OPTICAL  DATA  AND  METEOROLOGICAL  EVENTS 

Optical  and  meteorological  data  for  the  Relay  Mirror  Experiment  (RME)  were  col¬ 
lected  between  1985-1988  in  Kihei,  Maui,  Hawaii.  Reviewing  the  composite  Kihei  op¬ 
tical  dataset  reveals  the  Maui  site  as  "mediocre"  to  "poor".  While  upward  and  downward 
trends  occur,  the  predominant  transverse  coherence  length  values  generally  stay  within 
the  50  to  70  mm  range. 

The  lack  of  extremely  good  values  (e.g.,  200-400  mm)  is  tied  to  the  turbulent,  3-5° 
C  boundary  layer  inversion  2-3  km  above  the  Kihei  site.  Except  in  the  winter,  north- 
northeasterly  trades  usually  cap  the  boundary  layer.  The  interface  between  the  4  to  7 
m/s  trades  and  the  boundary  layer  provides  the  mechanical  turbulence  detrimental  to 
optimum  "seeing"  conditions  for  light  propagation.  The  concept  of  the  boundary  layer 
introducing  optical  turbulence  is  supported  by  Gossard  (1985).  Using  fast-response 
turbulence  sensors  traversing  a  300m  tower,  he  shows  a  turbulent  layer  of  approximately 
0.01  f?nr2li  to  exist  just  below  the  inversion.  This  implies  that  the  greatest  mixing  of  air 
masses  with  different  temperatures  (turbulence)  is  at  the  base  of  the  boundary  layer  in¬ 
version. 

Kihei  has  a  land/sea  breeze  even  during  the  "trade  wind"  season  (see  Chapter  2, 
Section  C).  This  daily  shift  of  circulation  generating  a  warm,  convective  air  mass  over 
Kihei  in  the  daytime  followed  by  a  cool,  relatively  dr}’  environment  at  night  renders  the 
mediocre  r0  measurements.  Katabatic  flow  off  of  the  3.1  km  high  Mt.  Haleakala  is  an¬ 
other  turbulence  producing  mechanism  observed  at  the  Kihei  site. 

Profiling  the  atmosphere  over  Kihei  during  summer  conditions,  one  expects  a 
land; sea  breeze  circulation,  a  boundary  layer  inversion,  trade  winds,  a  second  moist  layer 
aloft  and  a  moderate  jet  stream.  The  shear  surrounding  any  jet  is  a  source  of  turbulence 
exhibited  in  the  r0  magnitudes.  Shear  provides  the  mechanical  mixing  necessary  for  op¬ 
tical  turbulence,  however,  without  the  compliment  of  a  strong  temperature  gradient,  the 
effect  on  the  optical  regime  is  weakened.  This  is  best  exemplified  by  reviewing  the 
summer  trade  season  optical  measurements  collected  when  an  upper-level  jet  is  travers¬ 
ing  the  site  (weak  temperature  gradient).  The  0O  measurements,  which  are  sensitive  to 
turbulence  in  the  300  to  200  mb  jet  stream  level,  are  rarely  below  3  n  rad.  Note:  3  n  rad 
is  considered  poor  optical  conditions.  The  exception  to  the  rule  for  Maui  occurs  when 
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a  synoptic-scale  event  crosses  overhead,  tightening  the  temperature  gradient  associated 
with  the  jet  (see  Appendix  D).  A  stronger  example  is  examining  the  influence  of  a  polar 
jet.  The  tighter  temperature  gradient  and  subsequently  stronger  winds  associated  with 
the  polar  jet  are  able  to  drive  0O  measurements  down  to  1  to  2  n  rads,  a  very  turbulent 
environment  optically.  An  example  of  this  was  in  the  10  December  1985  Haleakala  data, 
where  0O  was  measured  as  1.6  n  rad. 

One  of  the  most  optically  active  Hawaiian  cases  occurs  in  April  1988,  a  winter  case. 
The  Kihei  site  experiences  the  passage  of  a  low  or  "eddy".  A  detailed  discussion  of  the 
optical  and  meteorological  conditions  throughout  that  experiment  is  found  in  Appendix 
D. 


B.  OPTICAL  VARIABLES  CALCULATED  FROM  RAWINSONDE  DATA 

During  four  of  six  RME  experiments  at  Kihei,  Maui  on-site  rawinsonde  launches 
were  released.  Two  sessions  employed  the  use  of  a  Vaisala  DigiCORA  MW  11 
Rawinsonde  System.  The  eight  April  1988  launches  and  ten  flights  recorded  in  July  1988 
began  with  a  two-second  sampling  rate  up  to  2  km,  followed  by  a  ten-second  sampling 
rate  for  the  remainder  of  the  flight.  All  18  vertical  profile  sets  were  utilized  for  testing 
four  distinct  Ot  calculation  methods.  The  final  rating  of  each  method  was  done  by  using 
the  Or  profile  to  calculate  the  two  optical  parameters,  r0  and  0O.  These  calculated  values 
were  compared  to  the  optical  measurements  acquired  simultaneously  with  the 
rawinsonde  launch. 

The  first  two  methods,  the  Simple  Potential  Temperature  Differences  (SPTD)  and 
the  Bulk  Gradient  (BG),  produced  similar  results.  Both  over  estimated  the  turbulence 
leading  to  an  under  estimation  of  r0  and  00  by  a  factor  of  10. 

Using  the  entire  profile,  the  1 1 1  and  13531  running  mean  techniques  produced  r0  and 
0O  values  of  the  same  order  magnitude  as  the  sensor  values.  The  r0  data  even  showed  a 
general  ability  to  match  the  relatively  large  and  small  magnitude  variations  observed  in 
the  sensor  data.  0o  calculated  from  the  running  mean  methods  were  generally  insensitive 
to  the  magnitude  changes  measured  by  the  isoplanometcr.  However,  the  calculated 
values  were  within  the  same  order  of  magnitude  as  the  sensor.  Overall,  the  1 353 1  run¬ 
ning  mean  provided  turbulence  parameters  at  least  16%  closer  to  the  sensor  values  than 
the  111  running  mean. 

Using  the  Richardson  number  to  isolate  the  turbulent  layer  from  the  non-turbulent 
brought  mixed  results.  Conditioning  the  Ot  profile  by  calculating  turbulence  for  layers 
with  Ri  <  0.25  generated  r0  two  to  three  times  larger  than  the  sensor.  The  BG  Case 
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provided  values  approximately  7%  closer  to  the  sensor.  The  running  mean  approaches 
were  within  8-71  mm  of  the  sensor.  00  calculations  for  SPTD,  BG  and  running  mean 
generated  results  such  that  about  40%  over  estimated  the  sensor's  0O ;  50%  under  esti¬ 
mated;  and,  10%  straddled  the  sensor  value. 

The  conditionally  unstable  results  (Ri  <  2)  for  the  SPTD  Case  under  estimated  the 
r0  sensor  by  17-80%. 

Poor  correlation  between  the  calculated  and  measured  optical  parameters  revealed 
the  existence  of  serious  problems.  Limitations  in  the  current  sensor/system  prevented 
successful  Cf  calculations.  The  limitations  were  broken  down  into  three  areas: 

1.  Excessive  Temperature  Probe  Response  Time, 

2.  Limited  Vertical  Resolution  (aliasing)  and 

3.  Limited  Temperature  Resolution  (Roundoff  Error). 

The  limited  temporal  resolution  of  the  DigiCORA  MW  1 1  rawinsonde  profiles  was 
confirmed  by  calculating  the  temperature  probe  time  constant  correction.  For  the 
2-sccond  data,  only  15%  of  the  actual  Cf  was  measured.  While  the  10-second  filter 
displayed  54%  of  the  actual  Cf,  this  was  deceptive.  The  distance  between  10-second 
samples  produced  aliasing  effects  when  the  mean  potential  temperature  gradient  was 
large.  This  inability  to  resolve  high  frequencies  artificially  enlarged  the  already  over  es¬ 
timated  turbulence  present.  In  addition,  this  large  40-45  m  separation  approached  the 
Kolmogorov  outerscale. 

The  2-second  data  was  not  as  susceptible  to  aliasing.  Instead,  the  0.1  K  temperature 
resolution  (roundoff  error)  introduced  quantization  noise.  Using  two  distinct  proce¬ 
dures,  the  standard  error  formula  and  a  vertical  Cf  quantization  noise  simulation  model, 
noise  corrections  of  6.2.vlO‘4A’J/?rJ/3  (2-second  data)  and  1.4.vlO_4A?/;rJ/3  (10-second  data) 
were  calculated.  Application  of  this  correction  confirmed  that  the  Cf  values  in  the  initial 
2  km  of  data  were  mostly  quantization  noise.  Because  most  optical  turbulence  is  ex¬ 
pected  near  the  surface,  this  realization  was  not  encouraging. 

Another  major  problem  with  calculating  Cf  from  rawinsonde  data  was  insufficient 
vertical  resolution  that  produced  a  fitting  problem  with  the  running  mean  case.  Just  as 
the  inability  to  resolve  higher  frequencies  (aliasing)  produced  an  over  estimation  of  Cf 
for  the  SPTD  Case,  this  effect  became  even  more  apparent  when  the  running  mean 
methods  encountered  a  3-5  6  C  inversion  with  a  10-second  sampling  rate.  This  enhanced 
the  existing  over  estimation  of  turbulence.  In  contrast,  turbulence  calculated  across  the 
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sharp  temperature  shift  (inversion)  captured  within  the  2-second  data  (1537  UTC  7  July 
1988)  did  not  appear  to  suffer  from  this  difficulty. 


VII.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


A.  LINKING  OPTICAL  DATA  AND  METEOROLOGICAL  EVENTS 

Quantizing  the  optical  turbulence  based  on  weather  conditions  is  a  non-trivial 
project.  As  simultaneous  optical  and  meteorological  observations  accumulate,  the  pos¬ 
sibility  of  plotting  "weather"  charts  of  iso  -  r0  and  iso  —  0O  lines  nears  reality.  Such  a 
task,  however,  requires  at  least  two  phases  of  investigation.  The  first  deals  with 
synoptic-scale  events.  Major  atmospheric  phenomena  afTecting  optical  propagation  in¬ 
clude:  the  passage  of  cold  and  warm  fronts,  the  polar  and  subtropical  jets,  high  and  low 
pressure  systems,  and  so  on.  A  dedicated  climatological  study  (spanning  many  years) 
of  optical  parameters  during  each  of  these  events,  with  measurements  taken  at  one  lo¬ 
cation,  would  provide  a  statistical  sample  from  which  ranges  of  values  could  be  deter¬ 
mined.  Being  realistic,  a  more  cost  effective  and  time  efficient  approach  for  acquiring  a 
(marginally)  valid  statistical  sample  is  to  design  a  set  of  experiments  with  multiple 
(minimum  3,  preferable  5)  r0  and  00  sensors.  The  geometric  shape  formed  by  4  of  the 
sensor  site  locations  would  be  a  square.  The  sides  of  this  square  would  be  along  the 
prevailing  winds  and  measure  a  kilometer  or  more  in  length.  The  fifth  sensor  would  be 
in  the  center  of  this  box.  Time  based  cross-sections  would  be  derived  from  three  of  the 
sensors.  Because  the  central  site  is  common  to  both  possible  cross-sections,  this  sensor 
would  be  the  most  critical. 

The  meteorological  sensors  needed  for  simultaneous  measurements  would  depend 
on  the  particular  event  being  studied.  For  example,  the  passage  of  a  cold  front  would 
require  surface  meteorological  stations,  rawinsondes  and  satellite  resources.  The  effect 
of  a  polar  jet  would  require  rawinsondes,  aircraft  and  satellite  data.  The  main  objective 
is  to  pinpoint  optical  measurements  (processed  in  real  time)  with  the  movement  of  the 
weather  phenomena.  Armed  with  at  least  three  "clean"  case  studies  of  a  specific  weather 
event,  r0  and  0a  ranges  for  various  sectors  of  the  specified  event  would  be  possible. 

The  second  area  of  investigation  involves  the  meso-  and  microscale  weather  activ¬ 
ities.  Though  weather  phenomena  on  these  scales  are  often  site  specific,  one  major 
source  of  optical  distortions  common  to  all  sites  is  the  boundary  layer.  The  influence 
of  the  boundary  layer  has  been  estimated  as  30-50%.  A  simultaneous  data  collection 
at  the  surface  and  above  the  boundary  layer  inversion  would  better  quantize  the  impor¬ 
tance  of  this  atmospheric  structure.  Because  Mt.  Haleakala  is  over  3  km  high,  such  a 
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project  could  theoretically  be  done  in  Maui.  However,  Haleakala's  Crater  experiences 
a  local  jet.  Therefore,  the  measurements  above  the  boundary'  layer  inversion  would  have 
to  be  sampled  sufficiently  far  from  this  accelerated  flow.  Another  possible  location  is 
along  the  central  California  coast  where  r0  values  have  often  been  recorded  to  be  more 
than  200  mm.  A  dedicated  study  contrasting  the  good  "seeing"  conditions  above  and 
below  the  boundary'  layer  interface  may  provide  a  better  resolution  for  this  question. 

Another  meteorological  phenomena  at  the  mesoscale  is  gravity  waves.  With  the 
current  method  of  data  acquisition  for  the  isoplanometer,  tracing  the  influence  of  gravity 
waves  in  real-time  is  possible.  This  ability  coupled  with  the  mobility  of  the 
isoplanometer  could  be  of  value  to  gravity  wave  studies. 

B.  CALCULATING  CJr  FROM  RAWINSONDE  DATA 

Another  approach  to  resolving  the  correlation  of  optical  turbulence  and 
meteorological  conditions  is  to  calculate  key  optical  parameters  (such  as  /•„  and  0O)  di¬ 
rectly  from  rawinsonde  data.  Currently,  one  of  the  main  obstacles  is  the  need  for  a  low 
cost  temperature  probe  and  operating  system  capable  of  acquiring  large  high  resolution 
data  bases.  An  adequate  single  probe  rawinsonde  temperature  sensor  would  require  a 
0.01-0.02  K  resolution,  100  msec  response  time  and  a  1  m  vertical  resolution.  These 
specifications  are  not  unreasonable  with  current  technology.  Some  method  of  detrend¬ 
ing  the  potential  temperature  profile  with  a  running  mean  still  will  be  necessary  until 
0. 1-0.2  m  vertical  resolution  is  possible. 

An  interesting  experiment  would  be  to  launch  a  rawinsonde  with  both  the  above 
recommended  single  temperature  probe  and  the  more  expensive  two  temperature  probe 
sensor.  The  latter  has  been  discussed  in  Olmstead  (198S).  A  comparison  cf  the  two  Cf 
profiles  generated  from  these  methods  would  help  to  calibrate  the  less  expensive  single 
probe  method. 

As  personal  computers  become  faster  and  high  data-acquisition  rate  rawinsonde 
packages  become  routine,  the  theoretical  concept  of  evaluating  and  predicting  the  op¬ 
tical  turbulence  over  any  site  in  real  time  will  evolve  and  become  a  reality. 


42 


APPENDIX  A.  CALCULATED/MEASURED  OPTICAL  PARAMETERS 

The  optical  variables  r0  (transverse  coherence  length)  and  6„  (isoplanatic  angle)  de¬ 
pend  on  an  integral  of  Q  along  the  vertical  path.  The  refractive  structure  parameter  Q 
is  proportional  to  the  temperature  structure  parameter  Q  which  can  be  computed  from 
the  rawinsonde  data.  Four  methods  for  computing  the  temperature  structure  parameter 
Ot  from  the  rawinsonde  data  include: 

1.  the  Simple  Potential  Temperature  Differences  (SPTD)  Equation, 

2.  detrending  the  potential  temperature  differences  with  a  111  Running  Mean  (111 
RM)  or  a  13531  Running  Mean  (13531  RM), 

3.  the  Bulk  Gradient  Equation  (BG)  and 

4.  conditional  data  sampling  based  on  the  Richardson  Number  (Ri<0.25  and 
Ri<2.0). 

Each  Or  method  was  tested  by  calculating  r0  and  0O  from  the  Q-  profile.  These  results 
were  ultimately  compared  to  the  optical  measurements  simultaneously  acquired  over  the 
same  site. 

Tabic  6  and  Table  7  correlates  the  flight  numbers  referenced  in  the  Appendix  A 
plots  with  their  rawinsonde  launch  time  (UTC). 


Table  6.  APRIL  1988  KIHEI,  MAUI  RAWINSONDE  LAUNCHES 


FLIGHT 

NO. 

TIME  (UTC) 

DATE  (DAY  MON  YR) 

1 

0610 

27  April  1988 

2 

1548 

27  April  1988 

0605 

28  April  1988 

4 

1549 

28  April  19S8 

5 

0605 

29  April  19S8 

6 

1558 

29  April  1988 

7 

0544 

30  April  1988 

8 

1434 

30  April  1988 
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Table  7.  JULY  1988  KIHEI,  MAUI  RAWINSONDE  LAUNCHES 


FLIGHT 

NO. 

TIME  (UTC) 

DATE  (DAY  MON  YR) 

1 

0703 

7  July  1988 

2 

1537 

7  July  1988 

3 

0536 

8  July  1988 

4 

1612 

8  July  1988 

5 

0605 

9  July  1988 

6 

1527 

9  July  1988 

7 

0556 

10  July  1988 

8 

1532 

10  July  198S 

9 

0550 

11  July  1988 

10 

0926 

11  July  1988 
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A.  R*  RESULTS 

Figure  9  through  Figure  11  display  the  r„  values  calculated  from  the  eight  April 
1988  and  ten  July  1988  Maui  rawinsonde  datasets.  An  estimate  of  the  r0  simultaneously 
measured  over  the  same  Kihei,  Maui  site  is  also  included. 

For  the  measured  r0  time  series  corresponding  to  these  profiles,  see  Nelson  (1988). 


45 


(mm)  r0  (mm) 


Transverse  Coherence  Length 


oo 
so 

70 
60 
50 
40 
30 
20 
10 
0 

0123450789 
Flight  Number  (April,  1988) 


S 


+ 

a 


A 


*****  SPTD 

A  AAA  A  BG 
Dooao  1 1  i  RM 
■H444  13631  RM 
Sensor 


A 


+ 

a 


A 

"T“ 


+ 

D 


A 

T“ 


■ 

+ 

D 


+ 

D 


+ 

O 


90  q 

SPTD 

80: 

aaaaaBG 

ooooo  1 1 1  RM 

70: 

+++++ 13631  RM 

■■■■■!>•  Sene  or 

00: 

■ 

50: 

m  m 

40^ 

4  *  4 

30  ^ 

D  *  i  *  s  *  0  ;  4 

20: 

■  D 

I0i 

********** 

0^ 

- 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

01  234507  80  10  11 


Flight  Number  (July,  1988) 


g.  9,  Measured  and  Calculated  r0  (All  Values  Included) 


46 


Transverse  Coherence  Length 


*****  SPTD 
*****  SO 

OOODO  1 1 1  RM 
+++++ 13631  RM 
Sensor 


Flight  Number  (April,  1908) 


( mm)  r„  ( mm) 


Transverse  Coherence  Length 

250 
200 
150 
100 
50 

0 

250 
200 
150 
100 
50 

0 

01  23450769  10  11 

Flight  Number  (July,  1980) 

Fig.  11.  Measured  and  Calculated  r0  (Ri  <  2.0) 


Flight  Number  (April,  1988) 


48 


B.  0,  RESULTS 

Figure  12  through  Figure  14  display  the  0O  values  calculated  from  the  eight  April 
1988  and  ten  July  1988  Maui  rawinsonde  datasets.  An  estimate  of  the  0O  simultaneously 
measured  over  the  same  Kihei,  Maui  site  is  also  included. 

For  the  measured  0O  time  series  corresponding  to  these  profiles,  see  Nelson  (1988). 
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Fig.  13.  Measured  and  Calculated  0  0  (Ri  <  0.25) 
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Fig.  14.  Measured  and  Calculated  0  0  (Ri  <  2.0) 
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APPENDIX  B.  MAUI,  HAWAII  -  AUGUST  1987 


A.  METEOROLOGICAL  DATA 

Appendix  B  includes  the  August  1987  vertical  profiles  of  atmospheric  data  acquired 
during  the  Relay  Mirror  Experiment  at  Kihei,  Maui.  The  data  acquisition  system  used 
was  a  VIZ  WO-8000  RP+  rawinsonde  system.  Vertical  resolution  is  approximately  35 
m.  Launch  site  was  at  20.45°  latitude,  156.27°  longitude.  Table  8  couples  the 
rawinsonde  flight  numbers  and  launch  times  for  all  the  August  19S7  Kihei  soundings. 


Table  8.  AUGUST  1987  KIHEI,  MAUI  RAWINSONDE  LAUNCHES 


FLIGHT 

NO. 

TIME  (UTC) 

DATE  (DAY  MON  YR) 

1 

0658 

18  August  1987 

2 

1440 

18  August  1987 

3 

0722 

19  August  1987 

4 

1103 

19  August  1987 

5 

1512 

19  August  1987 

6 

0607 

20  August  1987 

7 

1011 

20  August  1987 

8 

0650 

21  August  1987 

9 

1106 

21  August  1987 

10 

1559 

21  August  1987 

11 

0541 

22  August  1987 

12 

1053 

22  August  1987 
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August  07  -  Flight  1 


Mixing  Ratio  (g/kg) 


Fig.  15.  Kihei,  Maui  Rau insonde  Data:  0658  UTC  18  August  1987 


August  07  -  Flight  2 


Mixing  Ratio  (g/kg) 


Wind  Direction  (deg) 

200  300  400  600  600  700 


Fig.  16.  Kihei,  Maui  Rawinsonde  Data:  1440  UTC  18  August  1987 


August  07  -  Flight  4 


Mixing  Ratio  (g/kg) 


Fig.  18.  Kiliei,  Maui  Rauinsontle  Data:  1103  UTC  19  August  1987 
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Prearure-A Ititude  (km)  Pressure— Attitude  (km) 


August  07  -  Flight  5 


Mixing  Ratio  (g/kg) 


Potential  Temperature  (K) 


Wind  Direction  (deg) 


Kiliei,  Maui  Rauinsonde  Data:  1512  UTC  19  August  1987 


August  07  -  Flight  6 


Mixing  Ratio  (g/kg) 


Fig.  20.  Kiliei,  Maui  Rauinsonde  Data:  0607  UTC  20  August  1987 
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August  87  -  Flight  7 

Mixing  Ratio  (g/kg) 
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Kihei.  Maui  Rauinsonde  Data:  1011  UTC  20  August  1987 


August  07  -  Flight  0 


Mixing  Ratio  (g/kg) 


Fig.  22.  Kibei,  Maui  Rauittsonde  Data:  0650  UTC  21  August  1987 
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August  07  -  Flight  9 


Mixing  Ratio  (g/kg) 


Potential  Temperature  (K) 


Wind  Direction  (deg) 
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1  23456780  10  11 

Wind  Speed  (m/e) 


Fig.  23.  Kiliei,  Maui  Rawinsonde  Data:  1106  UTC  21  August  1987 


August  07  -  Flight  10 


Mixing  Ratio  ( g/kg ) 


Fig.  24.  Kihei,  Maui  Ran  insonde  Data:  1559  UTC  21  August  1987 
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August  87  -  Flight  11 


Mixing  Ratio  (g/kg) 


Fig.  25.  Kilici,  Maui  Ran  insonde  Data:  05-41  UTC  22  August  1987 
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August  07  -  Flight  12 


Mixing  Ratio  (g/kg) 


Potential  Temperature  (K) 


Wind  Direction  (deg) 


Kihei,  Maui  Rawinsonde  Data:  1053  UTC  22  August  1987 


APPENDIX  C.  MAUI,  HAWAII  -  DECEMBER  1987 


A.  METEOROLOGICAL  DATA 

Appendix  C  includes  vertical  profiles  from  the  December  1987  Kihei,  Maui 
rawinsonde  launches.  The  VIZ  WO-SOOO  RP  +  rawinsonde  system  used  to  acquire  this 
data  sampled  values  approximately  every'  35  m.  Launch  site  for  these  profiles  was  at 
20.54°  latitude,  156.27°  longitude.  Table  9  tabulates  all  the  flight  numbers  and  launch 
times  for  the  December  1987  Kihei  soundings. 


Table  9.  DECEMBER  1987  KIHEI,  MAUI  RAWINSONDE  LAUNCHES 


FLIGHT 

NO. 

TIME  (UTC) 

DATE  (DAY  MON  YR) 

1 

0950 

5  December  1987 

2 

1552 

5  December  1987 

0639 

6  December  1987 

4 

1131 

6  December  19S7 

5 

0517 

7  December  1987 

6 

No  Data 

- 

7 

1624 

7  December  1987 

8 

0537 

8  December  1987 

9 

1031 

8  December  1987 

10 

1557 

8  December  1987 

1) 

0345 

9  December  1987 

12 

095S 

9  December  19S7 
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Fig.  27.  Ki'ici,  Maui  Rauin  sonde  Data:  0950  UTC  5  December  1987 


Pressure— Altitude  (Jem.) 


December  07  -  Flight  2 


Fig.  28.  Kilici,  Maui  Rnuinsonde  Data:  1552  UTC  5  December  1987 


December  87  -  Flight  3 


Mixing  Ratio  (g/kg) 


Wind  Direction  (deg) 


Fig.  29.  Kiliei,  Maui  Ran  insonde  Data:  0639  UTC  6  December  1987 


Pressure-Altitude  (km)  Pressure-Altitude  (km) 


December  07  —  Flight  4 

Mixing  Ratio  (g/kg) 


Wind  Direction  (deg) 


Wind  Speed  (m/a) 


Fig.  30.  Kihei,  Maui  Rau  insonde  Data:  1131  UTC  6  Dcccml)cr  1987 


December  07  —  Flight  5 

Mixing  Ratio  ( g/kg ) 


Potential  Temperature  (K) 
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Fig.  31.  Kiliei,  Maui  Rauinsonde  Data:  0517  UTC  7  December  1987 


December  87  -  Flight  7 


Mixing  Ratio  ("'/kg) 
0123466709  10  11  12 


Wind  Direction  (deg) 


Fig.  32.  Kihei,  Maui  Rauin  sonde  Data;  1624  UTC  7  December  1987 
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December  07  -  Flight  8 


Mixing  Ratio  (g/kg) 


Wind  Direction  (deg) 


Fig.  33.  Kiliei,  Maui  Rauinsonde  Data:  0537  UTC  8  December  1987 
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Fig.  34.  Kihei,  Maui  Rauiusondc  Data:  1031  UTC  8  December  1987 


Kihei,  Maui  Ran  insonde  Data:  1557  UTC  8  December  1987 


December  07  —  Flight  11 


Mixing  Ratio  ( g/kg ) 


Wind  Direction  (deg) 


Kihei,  Maui  ft  an  in. sonde  Data:  0345  UTC  9  December  1907 


Pressure— Altitude  (km) 


December  07  -  Flight  12 


Mixing  Ratio  (g/kg) 


Wind  Direction  (deg) 


Fig.  37.  Kihei,  Maui  Rawinsonde  Data:  0958  UTC  9  December  1987 


APPENDIX  D.  MAUI,  HAWAII  -  APRIL  1988 
A.  METEOROLOGICAL  DATA 

Appendix  D  includes  the  April  19S8  vertical  profiles  of  atmospheric  data  acquired 
for  the  Relay  Mirror  Experiment  at  Kihei,  Maui.  The  data  acquisition  system  used  was 
a  Vaisala  DigiCORA  MW  11  rawinsonde  system.  Launch  site  was  at  20.75°  latitude, 
156.43°  longitude.  The  initial  2  km  were  sampled  every  two  seconds  (8-10  m).  Data 
acquisition  above  2  km  was  once  ever)'  ten  seconds  (40-45  m).  Flight  number  and 
launch  times  for  all  the  April  1988  soundings  are  provided  in  Table  10  . 


Table  10.  APRIL  1988  KIHEI,  MAUI  RAWINSONDE  LAUNCHES 


FLIGHT 

NO. 

TIME  (UTC) 

DATE  (DAY  MON  YR) 

1 

0610 

27  April  1988 

2 

1548 

27  April  1988 

3 

0605 

28  April  1988 

4 

1549 

28  April  19S8 

5 

0605 

29  April  1988 

6 

1558 

29  April  1988 

7 

0544 

30  April  1988 

8 

1434 

30  April  198S 
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Fig.  38.  Kihei,  Maui  Rauinsonde  Data:  0610  UTC  27  April  1988 


Flight  2 


Fig.  41.  Kihei,  Maui  Ran  insonde  Data:  1548  UTC  27  April  1988 
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April  00  -  Flight  2 

Mixing  Ratio  (g/kg) 


Kihei,  Maui  Ran  insonde  Data:  1548  UTC  27  April  1988 


April  BO  -  Flight  2 

Wind  Direction  (deg) 


Fig.  43.  Kilici,  Maui  Raumsoiide  Data:  1548  UTC  27  April  1988 
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44. 


Kihei,  Maui  Ran  in. sonde  Data:  0605  UTC  28  April  1988 


April  00  -  Flight  3 

Wind  Direction  (deg) 


Fig.  46.  Kiliei,  Maui  Rau  insonde  Data:  0605  UTC  28  April  1988 
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Fig.  *17.  Kiliei,  Maui  Raninsonde  Data:  1549  UTC  28 -April  1988 
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April  00  -  Flight  4 

M ixing  Ratio  (g/tcg) 


Fig.  48.  Kilici,  Maui  Rawinsonde  Data:  1549  UTC  28  April  1988 
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April  1U3  -  Flight  4 

Wind  Direction  (deg) 


Fig.  49.  Kilici,  Maui  RaWiusonde  Data:  1549  IJTC  28  April  1988 
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Kihei,  Maui  Rjwinsonde  Data:  0605  IJTC  29  April  1988 


April  08  -  Flight  5 
Wind  Direction  (deg) 


Wind  Speed  (m/s) 


Kihei,  Maui  Rawinsonde  Data:  0605  UTC  29  April  1988 


Pressure— Altitude  (km) 


April  00  -  Flight  6 

Mixing  Ratio  (g/kg) 


)  ~  |  h  i  n  i  m  rfrrrm  it  tjttttt  » i  npTrn  rrrrfn  m  i  TTTjrrrrn  i  ri^  TfYrn  m'|-  0 

200  250  300  350  400  450  500  550 

Potential  Temperature  (K) 


Kihei,  Maui  Rawinsonde  Data:  1558  UTC  29  April  1988 


Pressure -Altitude  (km) 


April  88  -  Flight  7 


Temperature  Sc  Dewpoint  (C) 


Fig.  56.  Kihci,  Maui  Rawinsonde  Data:  054-1  UTC  30  April  1988 
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April  00  -  Flight  7 

Mixing  Ratio  (g/kg) 


Potential  Temperature  (K) 


Kihei,  Maui  Ran  insonde  Data:  0544  UTC  30  April  I98fi 


April  00  -  Flight  0 

Mixing  Ratio  (g/kg) 


Tig.  60.  Kihei.  Maui  Rawinsonde  Data:  1434  UI'C  30  April  1988 
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B.  APRIL  1988,  EXPERIMENT  REPORT 

The  following  is  a  general  meteorological  summary  corresponding  to  the  optical 
data  measured  during  the  27  April  1988  (UTC)  through  01  May  1988  (UTC)  Kihei, 
Maui,  Hawaii  Experiment. 

1.  0600-1730  UTC  27  April  1988 

Extensive  cloud  cover  prohibited  optical  measurements  until  approximately 
1100  UTC.  The  0610  UTC  Kihei  rawinsonde  launch  includes  the  following  surface  pa¬ 
rameters: 

P  (Pressure)  =  1005.9  mb 
T  (Temperature)  =  22.6°  C 
RH  (Relative  Humidity)  =  71% 

Wind  Speed,  Wind  Direction  =  calm(<  2m/s),  N/A 
Stability  =  unstable  between  surface  and  0.2  km 

The  stratus  deck  restricting  optical  measurements  occurs  between  2.2-2.7  km,  and  ex¬ 
tends  horizontally  to  the  NE  3-5°  latitude  according  to  the  GOES-WEST  Infrared  (IR) 
Satellite  Image  taken  witnin  30  minutes  of  balloon  launch  time.  A  temperature  inversion 
of  5.1°  C  between  2.9  and  3.2  km  caps  the  cloud  layer,  and  is  coupled  with  light  winds 
( <  5  m  s)  from  the  NE. 

Winds  gradually  increase  with  height  until  9.9  km  (289.7  mb)  where  a  maximum 
wind  speed  of  22.0  m.'s  is  reported.  Wind  direction  remains  from  the  eastern  quadrants 
(0-180°  from  the  North)  up  to  13.2  km  (173.1  mb).  Within  the  next  97m  (174.4-171.1 
mb)  moderate  winds  become  calm  and  veer  201  degrees.  The  subsequent  increase  in 
wind  speed  peaks  at  13.7  m.'s  (143.2  mb)  before  it  gradually  decreases  to  winds  of  less 
than  10  m's.  Wind  direction  continues  to  be  from  the  NW  until  18.7  km  (well  after  the 
initial  signature  of  the  tropopause)  where  it  crosses  into  the  eastern  quadrants  for  the 
remainder  of  the  vertical  sampling.  The  rawinsonde  flight  is  terminated  after  90  minutes 
10  seconds  at  25.019  km  (24.9  mb). 

Comparing  the  0546  UTC  and  the  1246  UTC  IR  Satellite  Images,  the  dense 
stratus  deck  moves  from  NE  of  the  site  to  SW.  The  1246  UTC  and  1446  UTC  IR  Sat¬ 
ellite  Images  display  very  little  change  in  the  intermittent  stratus  deck,  though  r0  meas¬ 
urements  indicate  a  1  cm  increase  over  the  same  time  interval. 

Shortly  before  sunrise  (1548  UTC  27  April  1988)  a  second  rawinsonde  launch 

records  the  following  surface  parameters: 

P  (Pressure)  =  1005.3  mb 
T  (Temperature)  =  18.8°  C 
RH  (Relative  Humidity)  -  87% 
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Wind  Speed,  Wind  Direction  =  calm  ( <  2m/s),  N/A 
Stability  =  stable  conditions 

The  first  of  two  inversions  is  found  between  the  surface  (0.083  km)  and  0.257  km.  A 
second  slight  temperature  increase  of  0.4°  C  (0.365-0.430  km)  is  capped  by  a  small  hu¬ 
midity  increase  of  13%  between  0.459-1.012  km.  The  small  temperature  inversion  is 
supported  by  0.9  m/s  increase  in  wind  speed.  Wind  direction  remains  constant  at  28°  . 
Another  increase  in  moisture  content  is  displayed  between  1.270-1.737  km.  The  RH 
maximum  of  83%  is  reported  at  1.737  km.  Winds  sampled  thus  far  are  light  and  from 
the  NE,  NNE.  For  a  thin  layer,  1.782-2.275  km,  the  light  winds  back  to  northwesterlies. 
A  sharp  decrease  of  52%  RH  (2.075-2.713  km)  overlaps  the  variable  winds.  Embedded 
within  this  layer  is  a  1.5°  C  inversion  between  2.511-2.824  km.  Except  for  a  slight 
moisture  increase  between  7.754-S.824  km,  the  remainder  of  the  sounding  displays  an 
extremely  dry  atmosphere.  At  9.520  km,  winds  reach  a  peak  velocity  of  17.7  m.'s.  At  a 
slightly  higher  level,  10.558-12.396  km,  wind  speed  decreases  11.9  m/s  and  veers  163  de¬ 
grees.  The  veering  continues  another  54  degrees  over  the  next  0.943  km.  The  temper¬ 
ature  gradient  associated  with  the  217  degree  wind  shift  is  -20°  C/2.781  km  (2.61°  C/km 
less  than  the  drv-adiabatic  lapse  rate).  The  tropopausc  begins  around  17.3  km. 

2.  0600-1730  UTC  28  April  198S 

Kihei  rawinsonde  surface  measurements  recorded  at  0605  UTC  include: 

P  (Pressure)  =  1005.3  mb 
T  (Temperature)  -  ^2.6°  C 
RH  (Relative  Hunuuity)  =  72% 

Wind  Speed,  Wind  Direction  =  calm  ( <  2m/s),  N/A 
Stability  =  unstable  between  surface  and  0.2  km 

RH  reaches  a  maximum  of  91%  just  prior  to  the  first  of  two  inversions.  The  initial  in¬ 
version  is  between  2.254-2.430  km  (778.4-762.0  mb).  This  2.3°  C  increase  is  coupled  with 
an  8  degree  veering  of  winds,  and  no  change  in  wind  speed.  The  second  inversion  is 
fragmented  into  two  parts.  The  total  layer  displays  a  l.l&degee  C  increase  over  0.4  km 
(2.793-3.193  km).  Winds  associated  with  this  stable  layer  increase  by  1  m.'s  and  veer  15 
degrees.  As  before,  the  inversion  is  preceded  by  a  slight  increase  (+  18%)  in  RH. 

Between  5.627  km  and  6.004  km  light  winds  veer  132  degrees  NEE  to  SW.  The 
* 

temperature  gradient  associated  with  this  large  wind  shift  is  -3C/.377  km  (1.842  C/km 
less  than  the  dry-adiabaiic  lapse  rate). 
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Another  large  veering  in  wind  direction  is  between  8.951-10.830  km.  This  203 
degree  change  is  coupled  with  a  9.9  m/s  increase  in  wind  speed  and  a  temperature  gra¬ 
dient  of -13.2'1.S79  km  (2.775  C/km  less  than  the  dry-adiabatic  alpse  rate). 

The  jet  stream  occurs  between  1 1.561-14.431  km  and  flows  from  the  NN\V,NW. 
Peak  winds  of  42.1  m/s  at  315°  are  evidient  at  12.493  km.  The  temperature  gradient 
across  the  jet  is  -24.7°  C/2.870  km  (-8.606°  C/km).  The  Lihue  rawinsonde  data  recorded 
about  5.25  hours  after  the  local  Kihei  launch  displays  the  jet  stream  between 
12.354-14.580  km  (wind  direction  across  the  jet  backs  from  292°  to  279° ).  Temperature 
gradient  associated  with  Lihue's  jet  stream  measurements  is  approximately  -19.0° 
C/2.652  km  (7.164°  C/km).  Peak  wind  is  reported  as  42.8  m/s,  286°  at  13.487  km;  42.8 
ill's,  283°  at  13.940  km  (Note:  /•„  is  5  cm  during  both  sampling  intervals  0605  and  1115 
UTC.) 

At  1549  UTC  28  April  1988,  surface  measurements  arc: 

P  (Pressure)  =  1005.3  mb 
T  (Temperature)  =  17.9°  C 
RH  (Relative  Humidity)  =  90% 

Wind  Speed.  Wind  Direction  =  calm  ( <  2m,s),  N/A 
Stability  =  stable 

A  2.3°  C  inversion,  as  well  as  a  14%  drop  in  RH,  marks  the  first  84  m  of  the  vertical 
profiles.  Winds  remain  light  and  variable  for  the  initial  2  km.  A  second  minor  inversion 
is  between  2.368  and  2.S07  km  displaying  a  temperature  gradient  of  0.5°  C/0.438  km 
(1.134  °  C/km).  The  easterly  winds  associated  with  this  event  increase  2.6  m/s  and  veer 
6  degrees. 

A  general  wind  direction  shift  from  eastern  quadrants  to  western  quadrants 
takes  place  around  500  mb  (6.137  km).  Note:  500  mb  is  often  labelled  a  "level  of  non¬ 
divergence".  The  Level  of  Non-Divergence  is  defined  as  a  "mid-troposphc-ric  surface 
separating  major  regions  of  horizontal  convergence  and  divergent  associated  with  the 
typical  vertical  structure  of  the  migrator}’  cyclonic  scale  weather  systems"  (Huschke, 
1970).  In  the  upper  half  of  the  profile  (above  500mb),  RH  increases  at  a  rate  of  36% 
in  2.664  km.  This  secondary  peak  of  49%  RH  is  measured  at  12.005  km.  Nestled  above 
the  humidity  influx  is  the  jet  stream,  12.765-14.448  km.  Peak  winds  of  47.0  m/s  arc  re¬ 
corded  at  13. 319-13. 361  km.  The  temperature  gradient  across  the  jet  stream  is  -12.3° 
C/ 1.683  km  (-7.308°  C'km). 


3.  0600-1730  T,TC  29  April  1988 

Optical  data  gathering  during  this  session  was  severly  restricted  by  cloud  cover 

from  a  synoptical  scale  weather  system.  Unlike  previous  nights,  the  0605  UTC  sounding 

reveals  the  near  surface  winds  to  be  from  the  SW.  (ThU  implies  the  presence  of  an 

anabatic  flow  between  Kihei  and  Haleakala.)  The  rest  of  the  surface  parameters  are: 

P  (Pressure)  =  1006.5  mb 
T  (Temperature)  =  24.0°  C 
RH  (Relative  Humidity)  =  75% 

Wind  Speed,  Wind  Direction  =  calm  ( <  2m/s),  N/A 
Stability  =  unstable  between  surface  and  0.2  km 

The  vertical  humidity  profile  displys  two  peaks.  The  first,  a  gradual  13%/km 
increase  climaxing  at  84%  between  1.261-1.330  km.  A  minor  0.1°  C/0.008  m  inversion 
caps  this  first  humidity  maximum.  The  second  peak  follows  a  sharp  73%;km  gradient 
that  sends  the  humidity  soaring  to  75%  at  8.976  km.  Most  significant  to  this  sampling, 
however,  is  the  wind  profile.  As  stated  before,  near  surface  winds  (0.219-1.682  km)  are 
light  and  SW.  A  backing  into  the  eastern  quadrants  follows  and  persists  between 
1.759-5.104  km.  A  0.5°  C  inversion  interrupts  the  atmospheric  structure  between 
2.017-2.211  km.  Above  5.1  km,  winds  veer  into  the  western  quadrants  where  they  re¬ 
main  for  the  rest  of  the  v  ertical  sounding.  The  jet  stream  begins  at  13.046  km  and  ex¬ 
tends  to  14.902  km.  The  maximum  velocity  sampled  is  39.7  m/s  (13.721-13.773  km). 
The  temperature  gradient  accompanying  the  jet  is  -11.4°  CT.856  km  (-6.142°  C/km). 

Using  the  0546  UTC  IR  Satellite  Image,  the  site  appears  in  the  area  of  a  comma 
cloud  where  one  would  expect  overrunning  and  warm  advcction;  the  warm  front  would 
be  to  the  south.  At  1 146  UTC,  the  site  appears  to  straddle  the  cold  front.  The  Lihuc 
sounding  launched  at  1115  UTC,  records  a  surface  pressure  of  approximately  3  mb 
greater  than  Kihei,  after  interpolating  to  S3  m  (height  of  the  Kihei  site).  Since  Lihue  is 
approximately  2°  north  and  3°  west  of  Kihei,  this  post-cold  front  characteristic  is  ex¬ 
pected.  Associated  with  the  approximate  time  of  cold  front  crossing  (and  all  the  ex¬ 
pected  tightening  of  gradients),  the  r0  sensor  displays  a  20+  mm  drop  (!)  within  the  same 
time  interval.  At  1300  UTC,  r0  measurements  indicate  a  15  mm  plunge  in  magnitude. 
The  corresponding  satellite  image  implies  the  passing  of  additional  cloud  lines  associated 
with  smaller  fractions  of  the  larger  cold  front  during  that  time  period.  The  1446  UTC 
image  displays  the  passing  of  the  most  extensive  cloud  cover  associated  with  the  cold 
front  in  the  comma  cloud.  By  1546  UTC  (12  minutes  prior  to  Kihei  sonde  launch  6),  the 
site  is  boardering  the  area  that  is  normally  (in  more  northerly  latitudes)  associated  with 
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the  dry  tongue,  "a  protrusion  of  relatively  dry  air  into  a  region  of  higher  moisture  con¬ 
tent"  (Huschke,  1970).  In  this  sector,  one  expects  cold  advection  and  sinking  air. 

The  1558  UTC  Kihei  rawinsondc  launch  records  the  following  surface  parame¬ 
ters: 

P  (Pressure)  =  1006.6  mb 
T  (Temperature)  =  18.0°  C 
RH  (Relative  Humidity)  =  86% 

Wind  Speed,  Wind  Direction  =  calm  ( <  2m/s),  N/A 
Stability  =  stable 

A  small  inversion  of  2.8°  C  over  57  m  initials  the  profile.  Winds  remain  light  and  vari¬ 
able  throughout  the  first  5  km.  Humidity  measurements  indicate  a  small  peak  of  72% 
at  0.561  km  followed  by  two  minor  inversions.  A  second  peak  of  71%  at  1.715  km 
precedes  a  small  0.5°  C  increase  between  1.768-1.815  km.  A  third  humidity  maximum 
of  62%  is  recorded  at  3.159  km  with  yet  another  maximum  of  69%  at  7.658  km.  The 
highest  humidity  measurement  (apart  from  the  surface)  occurs  at  9.307  km  (73%),  with 
a  second  strong  influx  of  humidity  at  10.151  km  (72%).  This  moist  layer  immediately 
precedes  a  westerly  jet  stream  (10.964-13.037  km).  Unlike  the  former  soundings,  this 
jet's  maximum  wind  speed  is  38.6  m,s.  The  temperature  gradient  across  the  jet  stream 
is  -13.6°  C/2.138  km  (-6.361°  C/km).  Winds  remain  out  of  the  western  quadrants  well 
past  the  initial  signature  of  the  tropopause  (17  km).  Not  until  20  km,  do  the  winds 
suddenly  become  light  and  variable. 

4.  0500-1730  UTC  30  April  1988 

The  0544  UTC  Kihei  launch  surface  parameters  include: 

P  (Pressure)  =  1007.9  mb 
T  (Temperature)  =  23.3°  C 
RH  (Relative  Humidity)  =  79% 

Wind  Speed,  Wind  Direction  =  light,  W 
Stability  =  stable 

Surface  stability  is  defined  by  a  0.5°  C  inversion  in  the  first  6  m  of  the  profile.  An  in¬ 
termittent  layer  of  stratus  begins  at  2.082  km  (96%  RH).  A  second  stable  layer  marked 
by  a  1.3°  C  increase  in  temperature  over  0.379  km  (3.430°  C/km)  caps  the  moist  layer. 
By  2.687  km  the  winds  pick  up  and  veer  from  the  NW  to  S£.  The  second  of  three  peaks 
in  the  humidity  is  recorded  at  7.113  km  (84%  RH).  The  third  humidity  maximum  at 
8.982  km  (75%)  immediately  follows  a  layer  of  strong  winds  (33. S  m/s)  from  the  SW. 
A  weak  0.6°  C/0.101  km  (5.941°  C/km)  inversion  ties  this  layer's  structure  together.  A 
slow  westerly  jet  stream  (max  wind  speed  of  36.Sk  rn.'s  at  12.972-13.021  km)  occurs  be- 
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tween  12.716  km  and  13.228  km.  This  is  supported  by  a  -2.2°  C/0.512  km  (-4.297° 
C/km)  temperature  gradient.  Despite  a  decrease  in  wind  speed  with  height,  winds  remain 
W,S\V  up  through  18  km.  Around  the  height  of  the  tropopause  (20.5  km)  winds  back 
into  the  eastern  quadrants  where  they  remain  until  the  sampling  is  terminated  (23,975 
km). 

Consulting  the  1115  UTC  Lihue  sounding,  as  well  as  the  1146  UTC  satellite 
image,  the  overall  structure  of  the  comma  cloud  is  now  fragmented.  Low  level  clouds 
are  intermittently  dispersed  around  the  Hawaiian  Islands.  A  moist  layer  between  1-2.1 
km  is  implied  from  Lihue's  sounding.  Surface  winds  are  moderate  to  light  and  from  the 
east  until  about  6  km  where  they  veer  and  increase  in  magnitude.  The  maximum  velocity 
recorded,  however,  is  34.0  m.'s  from  the  west  (12.536  km).  This  is  hardly  the  strength 
of  the  system  just  monitored  less  than  48  hours  ago.  The  r0  measurements  reflect  this 
lack  of  gross  turbulence  with  an  increase  from  45  mm  to  nearly  70  mm  sampled  over  the 
same  5-6  hour  time  period. 

The  last  Kihei  rawinsonde  launch  of  the  experiment  (1434  UTC)  records  the 

following  surface  measurements: 

P  (Pressure)  =  1007.4  mb 
T  (Temperature)  -  18.0°  C 
RH  (Relative  Humidity)  =  100% 

Wind  Speed,  Wind  Direction  =  light,  NE 
Stability  =  stable 

The  vertical  profile  begins  with  a  3.4°  C  inversion  in  the  first  36  m.  Surface  humidity 
of  100%  (a  measurement  confirmed  by  condensation  on  the  telescope  lenses)  decreases 
20%  in  the  initial  0.312km.  The  light  NE  surface  winds  back  to  light  W,S\V  within 
0.265km  of  the  surface.  Humidity  increases  to  a  peak  of  83%  at  0.638  km,  dips  about 
5%,  then  increases  to  S7%  by  1.241  km.  A  0.5°  C  increase  occurs  between  1.324-1.440 
km  (4.310°  Okm),  followed  by  another  inversion  of  0.6°  C  between  2.315-2.464  km. 
Winds  display  an  increasing  trend  until  approximately  4.6  km  (SW).  At  this  same  level, 
humidity  plunges  to  values  of  less  than  10%.  The  westerly  jet  stream  is  nestled  between 
12.060-13.196  km  carrying  a  peak  velocity  of  37.9  m/s  betwen  12.444-12.521  km.  The 
tropopause  occurs  around  19.1  km. 

5.  0600-1730  UTC  01  May  1988 

Early  evening  r0  data  indicates  intermittent  cloud  cover.  The  0546  UTC  satellite 
image  displays  cloud  cover  as  a  local  event.  The  closest  cyclone  is  NW  of  the  site  by 
approximately  10°  latitude  (1100  km).  The  storm  system  that  traveled  through  the  area 
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earlier  in  the  week  has  dissolved  into  the  high  pressure  ridge  NE  of  Hawaii.  The  two 
available  soundings  (Hilo  and  Lihue)  launched  at  1100  UTC  and  1115  UTC,  respec¬ 
tively,  correspond  to  relatively  high  r0  measurements  (80  mm).  Both  rawinsonde 
soundings  indicate  an  absence  of  surface  inversion.  Dewpoint  depression  minimums  (no 
RH  is  recorded  in  these  two  rawinsonde  datasets)  occur  between  0.4  and  1.6  km  implying 
a  moist  layer  within  this  1.2  km  region.  Super-adiabatic  lavpse  rates  are  noted  between 
surface-1000  mb  (12.327°  C/km)  at  Lihue,  and  709-704  mb  (46.867°  C/km)  at  Hilo. 
Because  there  is  no  consistent  level  for  these  extremely  large  lapse  rates,  they  are  most 
likely  locally  induced.  Winds  below  4  km  are  generally  from  the  eastern  quadrants 
(Lihue's  easterly  winds  extend  to  5  km).  Above  this  level,  wind  speed  strengthens  and 
wind  direction  is  from  the  west.  Jet  stream  velocities  are  recorded  as  W,  NW  between 
11.792-13.922  km  over  Hilo,  and  W,SW  between  10.645-14.672  km  over  Lihue.  Maxi¬ 
mum  wind  speeds  of  42.8  m/s  at  11.785  km  (Hilo)  and  39.2  m/s  at  12.016  km  (Lihue) 
are  recorded.  Both  sites  measure  a  lapse  rate  of  -6.7°  C/km  across  the  jet.  The  two 
soundings  also  concur  that  the  tropopasue  caps  the  jet's  subsiding  winds  at  approxi¬ 
mately  17.3  km. 

r0  measurements  record  a  gradual  30  mm  decrease  between  1200-1300  UTC. 
The  only  available  resource  over  that  time  interval  is  satellite  images.  Apart  from  an 
apparent  increase  in  local  cloud  cover,  there  is  no  glaring  evidence  of  turbulent  sources. 
A  review  of  1988  May  2  reveals  the  development  of  another  comma  cloud  structure  west 
of  Hawaii.  Perhaps  the  sudden  drop  in  r0  (turbulence  increase)  is  a  results  of  the  over¬ 
running  (warm  advcction)  often  associated  with  the  NE  sector  of  this  synoptic  scale 
event. 
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APPENDIX  E.  MAUI,  HAWAII  -  JULY  1988 


A.  METEOROLOGICAL  DATA 

.  Appendix  E  includes  profiles  from  the  July  1988  Kihei,  Maui  rawinsonde  dataset. 
A  Vaisala  DigiCORA  MW  11  rawinsonde  system  was  used  as  the  data  acquisition  sys¬ 
tem.  Launch  site  for  this  session  was  located  at  20.75°  latitude,  156.43°  longitude.  The 
initial  2  km  were  sampled  every  two  seconds  (8-10  m).  Above  2  km,  the  sampling  rate 
was  once  every  ten  seconds  (40-45  m).  Flight  number  and  launch  times  are  provided  in 
Table  11. 


Table  11.  JULY  1988  KIHEI,  MAUI  RAWINSONDE  LAUNCHES 


FLIGHT 

NO. 

TIME  (UTC) 

DATE  (DAY  MON  YR) 

1 

0703 

7  July  I9SS 

2 

1537 

7  July  1988 

3 

0536 

8  July  198S 

4 

1612 

8  July  1988 

5 

0605 

9  July  1988 

6 

1527 

9  July  1988 

7 

0556 

10  July  19S8 

8 

1532 

10  July  1988 

9 

0550 

11  July  19S8 

10 

0926 

1 1  July  1988 
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Fig.  64.  Kihei,  Maui 


July  88  -  Flight  1 

Wind  Direction  (deg) 


sonde  Data:  0703  UTC  7  July  1988 
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Pressure— Altitude  (km) 


July  88  -  Flight  2 


Temperature  Sc  Dewpoint  (C) 


Fig.  65.  Kihei,  Maui  Ran  insonde  Data:  1537  UTC  7  July  1988 


July  08  -  Flight  2 

Mixing  Ratio  (g/kg) 


Fig,  66.  Kihei,  Maui  Raninsonde  Data:  1537  UTC  7  July  1988 


July  88  -  Flight  3 


Maui  Rawinsonde  Data:  0536  UTC  8  July  1988 
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July  80  -  Flight  3 

Mixing  Ratio  (g/kg) 


.  Kihei,  Maui  Rawinsonde  Data:  0536  UTC  8  July  1988 


July  08  -  Flight  3 


Wind  Direction  (deg) 


Fig.  70.  Kihei,  Maui  Rawinsonde  Data:  0536  UTC  8  July  1988 
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July  88  -  Flight  4 


Fig.  71.  Kihei,  Maui  Rauinsonde  Data:  1612  UTC  8  July  1988 
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July  88  -  Flight  4 

Mixing  Kitio  (g/kg) 


Fig.  72.  Kihei,  Maul  Rawinsonde  Data:  1612  UTC  8  July  1988 
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Fig.  73.  Kihei,  Maui  Ran  insonde  Data:  1612  UTC  8  July  1988 
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July  00  -  Flight  5 


Fig.  74.  Kihei,  Maui  Rawinsonde  Data:  0605  UTC  9  July  1988 
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July  88  -  Flight  5 

Mixing  Ratio  (g/kg) 


Fig.  75.  Kihei,  Maui  Ran  in. sonde  Data:  0605  UTC  9  July  1988 
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July  88  -  Flight  5 

Wind  Direction  (deg) 


Fig.  76.  Kihei,  Maui  Rauinsonde  Data:  0605  UTC  9  July  1988 


125 


July  88  -  Flight  6 


77.  Kihei,  Maui  Rawinsonde  Data:  1527  UTC  9  July  1988 


July  80  -  Flight  6 

Mixing  Ratio  (g/kg) 


78.  Kihei,  Maui  Rnuinsonde  Data:  1527  UTC  9  July  1988 


July  88  -  Flight  6 

Wind  Direction  (deg) 


Kihei,  Maui  Rauinsonde  Data:  1527  UTC  9  July  1988 
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July  88  -  Flight  7 

Mixing  Ratio  (g/kg) 
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Kiliei,  Maui  Rawinsomle  Data:  0556  U  I'C  10  July  1988 


July  88  -  Flight  7 


Fig.  82.  Kihei,  Maui  Rawinsonde  Data:  0556  UTC  10  July  1988 
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Pressure— Altitude  (km) 


July  88  -  Flight  8 


Temperature  Sc  Dewpoint  (C) 


§• 

Fig.  83.  Kihei,  Maui  Rawiusonde  Data:  1532  UTC  10  July  1988 
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July  88  -  Flight  8 

Mixing  Ratio  (g/tcg) 
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Kihei,  Maui  Rawinsonde  Data:  1532  UTC  10  July  1988 


Pressure  -A  Ititude  ( km ) 
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Fig.  85.  Kihei,  Maui  Rawinsonde  Data:  1532  UTC  10  July  1988 
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Pressure  —A  Ititude  ( km ) 


July  88  -  Flight  9 
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Fig.  87.  Kiliei,  Maui  Ranmsonde  Data:  0550  UTC  11  July  1988 


136 


-120  -100  -80  -60  -40  -20  0  20 

Temperature  &  Dewpoint  (C) 


rig.  89.  Kiliei,  Maui  Raivinsonde  Data:  0926  UTC  11  July  1988 
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July  08  -  Flight  10 

Mixing  Ratio  (g/kg) 


Fig.  90.  Kihei,  Maui  Rawinsonde  Data:  0926  UTC  11  July  1988 
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July  88  -  Flight  10 

Wind  Direction  (deg) 


0  50  100  150  200  250 


Kihei,  Maui  Rauinsonde  Data:  0926  UTC  1 1-  July  198.8 
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